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Lewis Research Center 
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SUMMARY 

This report describes OIGTEM, a digital computer program that simulates 
two-spool, two-stream turbofan engines. The turbofan engine model In OIGTEM 
contains steady-state performance maps for all of the components and has con- 
trol volumes where continuity and energy balances are maintained. Rotor 
dynamics and duct momentum dynamics are also Included. Altogether there are 
16 state variables and state equations. DIGTEM features a backward-difference 
Integration scheme for Integrating stiff systems. It "trims" the model equa- 
tions to match a prescribed design point by calculating correction coefficients 
that balance out the dynamic equations. It uses the same coefficients at off- 
design points and Iterates to a balanced engine condition. 

Transients can also be run. They are generated by defining controls as a 
function of time (open-loop control) In a user-written subroutine (TMftSP). 
DIGTEM has run on the IBM 370/3033 computer using Implicit integration .with 
time steps ranging from 1.0 msec to 1.0 sec. 

DIGTEM Is generalized In the aerothermodynamlc treatment of components. 
This feature along with DIGTEM' s "trimming" calculations at a design point 
makes It a very useful tool for developing models of specific engines having 
the same two-spool, two-stream configuration. Also subsets of the turbofan 
engine configuration such as a turbojet or a turboshaft can be simulated with 
minor modifications to the Fortran coding. With extensive modifications to 
the coding, arbitrary configurations can be modeled. 

Included In this report Is complete documentation of DIGTEM. Input 
requirements, flow charts, modeling equations, and a test case are given along 
with a listing of the user-written subroutine TMRSP. Finally the use of DIGTEM 
to generate models for engines whose configurations are subsets of the general- 
ized turbofan engine configuration Is described. 


INTRODUCTION 

The development of aircraft propulsion systems depends, to a great extent, 
on being able to predict the performance of the propulsion system and Its 
associated controls. Computer simulations provide the means for analyzing the 
behavior and Interactions of these complex systems prior to the building and 
testing of expensive hardware. Simulations can also serve as aids In under- 
standing and solving problems that arise after the propulsion system Is 
developed. 


Computer simulations can be either generalized or specific to a particular 
propulsion system. Generalized simulations are desirable In that they allow 
for paper studies of many different engine configurations. Many generalized 
digital engine simulations exist today. Most of them are limited to steady- 
state performance calculations for a fixed number of engine configurations 
(refs. 1 to 4). But, since they are generalized, the user need only specify 
which of the configurations Is to be analyzed and supply the correct Input 
data. One generalized code, NNEP (ref. 5), lets the user build up arbitrary 
configurations through Input definitions. Another generalized code, DYNGEN 
(ref. 6), has transient capability but Is limited to the fixed engine con- 
figurations of references 3 and 4 (GENENG and GENENG II). All of the general- 
ized codes described have various limitations: they are limited to steady- 
state calculations, or they have many but fixed engine configurations. Some 
(DYNGEN and GENENG) are difficult to change, and none can scale Its model 
equations to reflect real engine data. Thus there Is a need for a computer 
code that can do both steady-state and dynamic calculations, Is flexible for 
modeling various engine configurations, and can also be easily adapted to model 
real engine data. 

Such a generalized dynamic engine program has been developed for the 
hybrid computer. That program Is called HYOES (ref. 7) and can handle the 
same fixed engine configurations as OYNGEN. However, by utilizing the capa- 
bilities of both the analog and digital computers, HYDES Is able to provide 
Improved engine model fidelity and an Interactive user environment. 

Even when using the HYDES program, the development of hybrid simulations 
Is time consuming and requires experience In dynamic system modeling, hybrid 
computer programming, and hybrid computer operations. To simplify this devel- 
opment process, a systematic, computer-aided approach for generating hybrid 
computer simulations of a particular class of engine (l.e., two-spool, two- 
stream turbofan) has ueen developed (refs. 8 and 9). This approach features 
more generalized aerothermodynamlc models of engine components and automated 
calculation of scale factors and simulation coefficients. Also a specified 
operating point, designated as the design point, Is used to scale the component 
maps and to determine correction coefficients that will balance the dynamic 
equations at the design point. This assures good steady-state accuracy at the 
design point. 

Despite the advantages of hybrid simulations they are generally not port- 
able or easily modified. Thus a digital computer model possessing the capa- 
bilities of the hybrid model presented In references 8 and 9 Is desirable. 

Such a model has been developed and Is the subject of this report. The digital 
portion of the hybrid model was retained and the dynamic equations that were on 
the analog were added to the digital code. The model was also unsealed to make 
It easier to modify or to Integrate with controls. A. numerical Integration .. 
scheme was added to provide dynamic capability to the digital program. The 
Integration technique Is Implicit and Is well suited for Integrating “stiff 11 
systems such as the turbofan engine model. The resultant digital computer code 
Is called DIGTEM for DIGItal Turbofan Engine Model. 

DIGTEM Is generalized In a different sense than DYNGEN. DIGTEM, although 
having only one engine configuration In the code, Is written In modular form to 
permit variations of the engine configuration (l.e., turbojets and turboshafts) 
to be simulated. This provides more flexibility (at the cost of recoding the 
Fortran) than DYNGEN, which Is limited to a fixed set of configurations and 
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Si 1 ?!*!! t( > cha ? 9e i Both OYNGEN and DIGTEM do component map scaling 

t? o r^SfJ?f?^JfJ a J ? de51 ?!! p 2 1nt \ However » 0I6TEM also calculates correct 
tlon coefficients to balance the dynamic equations so that a steady-state 

rJlfSS® des1 9" Point Is generated. The same values of the correction 

coefficients are used at off-design points. If the coefficients do not balance 
hJ? a 5 y Tl C f* uat1 °ns ft the operating points, DIGTEM Iterates to a new 
balanced eng ne condition. OIGTEM's flexibility should allow It to be a useful 
tool for engine dynamics studies and controls analysis. 

ranoiSJ 6 fr!Im h n S i b !Sfl3 r !! n ? ft n the IBW f 70 ^ 3033 Mainframe computer with time steps 
[fIJ Sd 7 sec t0 J'° s ! c ‘ s1nce 0I6TEM utilizes an Implicit Integra- 

t on scheme, the larger time ster-, can be used to generate fast, stable, tran- 

llfii s ? lut 1f n s. However, If the time step Is too large relative to the 
smallest engine time constant, there will be a loss In dynamic accuracy. 

mentJhfL^fSI^r 0 ^®^ do ;j m * ntat1 °" of the DIGTEM program-. Input requlre- 
TnJi, S ,JiJi 2k 2 h i r J s ’ and modeling equations are provided. A test case Is 
Included that makes use of a user-written, open-loop control subroutine, TMRSP. 

i]+hVn«laTr P m*? use ff 15 provided as a section of this report. Check 

with COSMIC, University of Georgia, Athens, Ga. 30602, for the availability of 
this program. Finally the use of OIGTEM to model turbojet and turboshaft Y 
engines is described* 


MODEL DESCRIPTION 

The engine model supplied with DIGTEM represents a two -spool, two-stream 
augmented turbofan engine. Figure 1 shows a schematic representation of that 

I fl w e 4 I 1 ®! 1S 2 Sed t0 supply airflow to the fan. Air leaving the 

InJ+hL 222c? ted d nt ° J W ° strea ? s ‘ one Passing through the engine core and 
another passing through an annular bypass duct. The fan Is driven by a low- 

hw e f S M2h tUrb1ne ‘ I h \f ore airflow passes through a compressor that Is driven 
by a high-pressure turbine. Both the fan and compressor are assumed to have 

a?r^t^?+2H e «+ y +h 0r be ! ter stab11 )J y at low speeds. Engine airflow bleeds 
are extracted at the compressor exit (station 3) and used for turbine coollna 

hn!^H r + tUrnS H t0 tbe , cyc1e >‘ Fuel flow Is Injected In the main combustor and 9 

^rM^ t 2 ft mK? dUC ? h0t 9aS fo I dr1v1ng the turbines. The engine core and bypass 
streams combine In an augmentor duct, where the flows are assumed to be 

thoroughly mixed. Additional fuel Is added to further Increase the aas tem- 
perature (and thus thrust). The augmentor flow Is discharged through a vari- 
able convergent-divergent nozzle. The nozzle throat area (station 8) and 

m?J?m5L n ^ Z V r ? d (stat1 °? E ) are varied to maintain engine airflow and to 
minimize drag during augmentor operation* 


Figure 2 contains a computational flow diagram of the engine model All 
ij*"’? *"> f fined in appeal, A. The analytic^ model 1n"Sd« SSniiartate 
maps to model the steady-state performance of the engine's rotating components. 
Fuld momentum In the bypass duct and the augmentor, mass and energy storage 
? control volumes, and rotor Inertias are Included In the model to provide 
transient capability. The complete engine model Is presented In appendix B. 

nm»ilr h u! n ^2 rat1 2! tecl »n1que used In DIGTEM Is a backward-difference 
(Implicit) Integration scheme that Is well suited for Integrating "stiff 

systems. A typical engine model will have time constants that differ by three 
or four orders of magnitude. This requires the use of very small time steps 
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It?" U Th2 B bS?k!2^‘d!f^™^ i?2 P I ,c1t> ""^ration schemes ta Insure stabtl- 

•nd W, 3i the six gas temperatures T s , T 4 , T 4 ., t T s , t£ Si "fjjj Si Si 
two mess flow rates_ w 13 and i 6 . The ordering of the state variables In 

l\rrv e ::° r VS ,S ^ rtant ' Tha Of VS and the corresponding 

follows? 1 " Va VG vecto VD0T are defined In terms of computer variables as 


VS(1) - XNL 
VS(2) . XNH 
VS(3) . W3 
VS( 4) - T3 
VS(5) - W4 
VS(6) - T4 
VS(7) . W41 
VS<8) . T41 
VS(9) « W6 
VS(10) . T6 
VS(11) » W7 
VS(12) = T7 
VS(13) « WA13 
VS(14) - WG6 
VS( 15) a W13 
VS(16) = T13 


VD0T(1) . OXNL 
VOOT(2) - OXNH 
VDOT (3) - DW3 
VDQT(4) » 0T3 
VDOT (5) > DW4 
VD0T(6) . DT4 
VDOT(7) - DW41 
VD0T(8) . DT41 
VD0T(9) a 0W6 
VDOT(IO) - DT6 
VDOT ( 77 ) „ DW7 
VDOT(12) * 0T7 
VDOT (13) « 0WA13 
VDOT(14) * DWG6 
VD0T(15) « 0W13 
V00T(16) » 0T13 


and tIS e cSr r e e n r o??ir?a??S?es M a r J? 6 l? S t^ end.^Thls^facilit^t ? u Sl! var,ables 


jte* though the Integration scheme featured In DIGTEM Is a hariruarH 

1 tel'Sjtlin? 'TotVllliril 'jjesretlon Shemi (Euler) 

OIGTEh Is described “ ttTSrtL Tmmm. ° pt ’ 0 " s 


USERS MANUAL 
Simulation flow Diagram 

flow chart for thl’ Ml^X^OlGTEIl! 8 lim/oiOTEiTSrltM out a^headlng 

4 


Identifying the type of engine being simulated. User data are then read In to 
define the Integration time step, printout Interval, operating point, and tran- 
sient duration. Next, INDATA Is called to read in component maps abd steady- 
state operating-point data. Both nonaugmented (dry) and augmented (wet) oper- 
ating points may be Input. By definition, the first dry point and the first 
wet point are design points. Once the desired design point or points Is speci- 
fied, DSGNPT Is called to calculate the scaling coefficients for the specified 
dry and wet design points. Next the engine parameters are calculated In ENGINE 
and vectors are set up for the Integration routines. Then BDINTG or FDINTG Is 
called to generate transient results depending on which Integration method Is 
desired. Once the transient Is completed, the simulation Is stopped. 

Flow charts for the subroutines are shown In appendix D. The following 
list defines the functions of the various subroutines: 

DIGTEM The main program for the simulation Is used to control the simulation. 
BDINTG The BDINTG subroutine performs Implicit Integration of the dynamic 
equations In DIGTEM. This subroutine Is discussed In detail In 
appendix C. 

BDPRNT The BDPRNT subroutine prints out either a short or a detailed output 
when backward difference Is used. 

DCTINT. The OCTINT subroutine calculates the derivative of the duct flow and 
performs a forward-difference Integration If desired. 

DMINV The DMINV subroutine performs a double-precision matrix Inversion of 
the Jacobian error matrix. 

DSGNPT The DSGNPT subroutine Is used to calculate correction coefficients from 
design-point data. At the dry and wet design points the scaling 
coefficients are calculated from Input values of pressure, tempera- 
ture, flow, etc. The correction coefficients are used to compensate 
for small modeling errors (e.g., map Interpolation errors or mis- 
matched component models) and to give zero derivatives at the design 
points. Additional coefficients are calculated at the wet design 
points so that a balanced condition exists In the augmentor at the 
wet design (maximum thrust) point. This subroutine Is discussed In 
more detail In appendix B. 

DUCT The DUCT subroutine calculates the duct Integration constants and 
losses. 

ENGINE The ENGINE subroutine solves the turbofan engine model by using the 
correction coefficients from DSGNPT and by calling the 14 engine 
subroutines In order. This routine Is called by DIGTEM to calculate 
Initial conditions when forward-difference Integration Is specified. 
ENGBD The ENGBD subroutine performs the same function as ENGINE but Is used 
when backward-difference Integration Is specified. 

ENG1 The ENG1 subroutine calculates fan performance. 

ENG2 The ENG2 subroutine calculates compressor performance. 

ENG3 The ENG3 subroutine calculates continuity and energy balances In the 

fan mixing volume and performs Integration If not In steady state. 
ENG4 The ENG4 subroutine calculates bleed flows. 

ENG5 The ENG5 subroutine calculates continuity and energy balances In the 

compressor. mixing volume and performs Integration If not In steady 
state. 

ENG6 The ENG6 subroutine calculates the high-pressure-turbine performance. 

ENG7 The ENG7 subroutine calculates the continuity and energy balances In 

the combustor and performs Integration If not In steady state. 

ENG8 The ENG8 subroutine calculates the low-pressure-turbine performance. 
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ENG9 The ENG9 subroutine calculates the continuity and energy balances In 
the high-pressure-turbine mixing volume and performs Integration If 
not In steady state. 

ENG1G The ENG1Q subroutine performs continuity and energy balances In the 

low-pressure-turbine mixing volume and performs Integration If not In 
steady state. 

ENG11 The ENG11 subroutine calculates nozzle performance. 

ENG12 The ENG12 subroutine calculates continuity and energy balances In the 

augmentor mixing volume and performs Integration If not In steady 
s tdto * 

ENG13 The ENG13 subroutine calculates the low- and high-spool derivatives and 
performs Integration If not In steady state. 

ENG14 The ENG14 subroutine calculates duct parameters and performs Integra- 
tion If not In steady state. 

ERROR The ERROR subroutine calculates the error vector for Implicit 
Integration. 

FOINTG The FDINTG subroutine performs forward-difference Integration of the 
dynamic equations. 

FLCONO The FLCOND subroutine calculates ambient pressure and fan Inlet total 
pressure and temperature from specified values of altitude, Mach 
number, and sea-level ambient temperature. 

FOOR The FOOR subroutine Indicates when data Input to a function Is out of 
ra^ge of the function data. It also Indicates If the Implicit 
Integration routine Is generating a Jacobian matrix. If MATRIX « 0, 
the routine Is not -generating a new matrix and the data should be 
checked. 

FUN! The FUN1 subroutine does a single-value Interpolation. 

FUN1L The FUN1L subroutine Is called following the cill to FUN1 when the same 

pair of breakpoints can be used to compute a second function value. 

GUESS The GUESS subroutine updaiftS..tlie guess vector for the Implicit Integra- 

tion scheme. 

INOATA The INDATA subroutine Initializes and reads In map data. 

MAP The MAP subroutine does a double-value Interpolation. 

MAPI The MAPI subroutine Is called following a call to MAP when the same 
four breakpoints can be used to- compute a second function value. 

MOOR The MOOR subroutine Indicates when data Input to a map are out of range 
of the map data. It also Indicates If the Implicit Integration 
routine Is generating a Jacobian matrix. If MATRIX * 0, the routine 
Is not generating a new matrix and the data should be checked. 

NOZZL The NOZZL subroutine calculates nozzle performance for both subsonic 
and supersonic flow conditions. 

PROCOM The PROCOM subroutine calculates the values of JP-4/alr thermodynamic 
properties based on supplied values of temperature and fuel-air 
ratio. The thermodynamic properties are the specific heats, the 
specific heat ratio, and the specific enthalpy. 

SPLINT The SPLINT subroutine calculates the spool speed derivative and 
performs forward-difference Integration If specified. 

SPOOL The SPOOL subroutine calculates the spool Integration constant from 
the moment of Inertia. 

TPRINT The TPRINT subroutine prints out short or detailed output when 
forward-difference Integration Is used. 

TRAT The TRAT subroutine calculates the Isentroplc temperature rise param- 
eter based on specified values of pressure ratio and specific heat 
ratio. 
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VQLINT 

VOLUME 


The VOUNT subroutine performs continuity and energy Integrations for 
forward-difference Integration and forms derivatives for the 
backward-difference Integration. 

Th the°1deal -Sas°1ai?° cakulates contro1 volul » stor<M ' •>» using 


A final subroutine TMRSP Is user written, 
function of time for transient operation. 


It defines open-loop controls as a 


program sotup 

Data are Input to DIGTEM via three methods. First, all comoonent man Hatn 
and operating-point data are specified In an Input data set* second P data 
1*^ ° H " routine selection. Integration time step, printout options, and 

MJ!ili? t aJ! t !^ ar !«! p ! c J f1ed 1n the ma1n r °vt1ne DIGTEM; and third, open-loop 
controls are specified In a user-written subroutine TMRSP. " 

fion S^ The 1n PUt data set supplied with DIGTEM Is shown In 

INDATA 4, ThI h f?r«+ d i?n are r ! a ? n by ! he ma1n Program DIGTEM and by subroutine 
INDA7A. The first line contains constants for the fraction of turbine coolino 

r ]rJL t J at + perfor ^ work for the h1 9 h “ a " d low-pressure turblnes^respectlvel^ 

Is (5F?2 t 5? t C ?he ? BLWHT and K C» WLT and the Input format' 

to dry design ^rX-pSlnt^ues *" C °"" > °" ent that ar * 

a ,en^:? r d e ,pSs^:‘ e :? s t;: d e £ 

Sn?T : b ? ?: e h r c °rr U 2° 

format They are Th f1fSt 11 f th data conta1ns f1ve numbers In (513) 


MAPNO NCV NPT NFCT NCOM 

SJS r ? M J!! N0 1s a map number to be used In the Za • f« (X,Y) function call* 

NCV is the number of curves Y on the map; NPT Is the number of oolnts tt 7 \ 
on each curve; NFCT Is the number of common functions l\ of the same * 
Independent variables; and NCOM Is the switch to Indicate that the* I 

point values can be used for all of the MCV cufSe (A !erS Indicates that 

the X values are different for each curve.) indicates that--- 

map statements 1 ^! t0 5\ used ^reading the remaining 

map statements, a (BX,7(4A2)) format Is used to specify the formats of m 

the X values, (2) the Y values, and (3) the z, values The reLlni™ 

JhS e ^l n + the data * u et S of,ta1n * 1n order ) the Y values, the X values for 
the first curve, the Z^ values for the first curve of each function the 

i~«s; 3 nxn m 

onc”l«nedlately a fol lowing' the < ? C °™Tu« . th °” * Val “ eS " Md be ° n1 » , '» lut 

uaMa 5?^ returning to figure 4, the first Input data set Is for the fan 

?ec tSd firSrtlL'IIr.; J5’ S ? ,ves the ad J ust "»"* to tde value oJ nor- 
recieo ran airflow due to off-schedule geometry* The efferf< aro mn^ni ^ 

the a iipSts? aP Thlt h 1s? n var1able ‘ 9eometr y Position and fan corrected speed as 
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AW,) 


(fW, 




The first line of data Is 


) 


0) 


7 ^ 1 0 (understood) 


EiS2 a«S*S 

values (normalized speeds), lines 5 and 6 2™ L nos * and 4 the v 

(the FVQP values - net norm zed! aJd ?w ? h ! * V,1 “°J f « r curve 1 
curve 1 (the flow shifts* /If?? a !25 ;l nes 7 u? nd 8 are the 2 values for 
for curve 2 , and lines 11 and 1? apf^+hn 1 ? are values (FVQP positions) 

12 more curves of flow shift as a functlon 0 of e rvrB d n 3 How shifts. Data for 
that the twelfth value of corrected Sid i! i oS a "?J ,p# J d fo " ow ’ Note 
ms curves (lines 45 to 52, pass throSgh the J«?g°n pl?"%r! ft.TITo)!' 

This 2 l gives $ the°sh1 f M ^corrected SlTf f^T U J] e - (l r etr » effects 
seometry. As Is the case of the fan ?he off -“hedule compressor 
actual variable-geometry position and correct peed” correlate with 


AW 


2.2 
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The first line of data Is 


(cvsp. 


( 2 ) 


14 


11 


1 


0 (understood) 


map number!” 11 The nexM 1 ne°1 s^aa 1 n 6 tho°f n t I y + da f a ‘ The 2 corres Ponds to the 
and 4 are the normalized JptN? ?n« s f °M ead1n 9 the data. Lines 3 

Thlrteen^moJe corrected^^eeMoDo)! ? 

are the normalized fan speeds V e^. Unes 5 and 5 are the fan duct 

po1nt U yalIe! 1 °L1nes 3 7 l> and 8 B def1n2 P the V ?orrected fan"?) 1 ” 11 *° *5® de5 ’ 9 "- 
corrected speed (0.3000): orrected fan flow curve for the first 


(wj 


fan,M 


(hi ! l \ 

X 2 ' «l /2 J 


(3) 


with fan^rfJrmSnce^ maps associated 


^ fan^OD 


f, 


(hi n l \ 

1 V* ’ °vj 


( 4 ) 
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the fan tip region pressure ratio. 



ORIGINAL PAGE m 
OF POOR QUALITY 

( 5 ) 


and finally the fan hub region efficiency 


"fan, 10 * 



( 6 ) 


The next set of data Is for the baseline (scheduled CV6P) compressor per- 
formance. Herv NFCT ■ 2. Thus there are two maps with the same Input vari- 
ables, namely pressure ratio and speed. Lines 3 and 4 are the normalized 

/ I /2 

2 * Lines S and 6 are the normalized compressor 

pressure ratios P 3 /P 2.2 for the first speed (0.700). Lines 7 and * contain 
the corresponding compressor corrected flows for the first speed: 



(7) 


The second map. for the compressor Is the compressor efficiency: 



( 8 ) 


The next set of data Is for the high-pressure-turbine performance. Here 
NFCT ■ 2. Thus there are two maps. In this case, NCOM * 1. Hence all curves 
for both maps are defined by the same pressure ratio breakpoints. Each map 
has eight curves with nine points per curve. Line 3 contains, the eight 

normalized speed values N^ Line 4 contains nine normalized pressure 

ratio values P3/P9. ?• Line 5 defines the-normallzed flow parameter curve 
for the first speed (0.7129): 



Line 6 defines the first curve of the map (l.e., the turbine enthalpy drop 
parameter) : 



(10) 


Lines 7 and 8 define the same two maps at the second speed (0.7548). Note that 
the same pressure ratio breakpoints are assumed for all speeds. 
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The final set of component map data Is for the low-pressure-turbine per- 
formance, The data are organized In exactly the same way as those for the 
high-pressure turbine, Lino 5 Is the normalized flow parameter for the first 
speed (0.4Q76) : 

N, 


<v 

p LT 


f 


16 




( 11 ) 


and line 6 Is the normalized turbine enthalpy drop parameter: 



After the component map data are read In, a blank line must bo Inserted. 
The next set of Input data are the bias values for the fan and compressor 
variable geometry. These values are subtracted from the actual variable- 
geometry position values In the simulation so that Inputs to the flow shift 
maps are always positive (a requirement of the Interpolation routine). The 
biases are 0.000 and 4.000 for BSFVGP and BSCV6P, respectively, and the format 
Is (2F10.0) . The next two numbers are the number of dry and wet operating 
points ( NDRY and NAU6, respectively) to be Input by the user. The format Is 
(IX, 2(12, 2X)). 

The following number Is a label for designating the operating point 
(POINT). The format Is (9X,I3). The first operating point Is assumed to be 
the dry design point. The operating-point data are organized as follows: 


P0 

P2 

PI 3 

P22 

P3 

P4 

P41 

P5 

P6 

P7 

TAM 

T2 

T13 

T22 

T3 

T4 

T41 

T6 

T7 

WA2 

WA13 

WA22 

WA3 

WG4 

WG41 

WG6 

WG7 

DH4 

0H41 

ETAB 

ETAAB 

FN 

XNL 

XNH 

WF4 

WF7 

A8 

AE 

ALT 

XMN 

CDN 

CVN 

FVGP 

CVGP 

FG 

The format Is (5F12.5). 

The next 

set of 

variables read 

In are 

volumes, reactances, and 

rotor moments of 

Inertia: 


VI 3 

V 3 

V 4 

V41 

V6 

AQL13 

AQL6 

XIH 

XIL 



the physical 


V7 


The format Is (6F12.5). The final values read In for the operating point are 
the fan tip region, fan hub region, and compressor efficiencies: 

ETAOf ETAIF ETAHC 


The format Is (3F12.5). 

The following sets of data contain additional dry-operating-point data. 
Once all the dry operating points are read In, the augmented operating points 
are read In. Again, the first wet operating point Is assumed to be the wet 
design point. 


10 


( rSI S } ^ r! fi ?£g c 1 M IS* ?i! ta f ? r “ ,B I E " *«■»»«*, 

data. are listed in I nqppr ri«S^l! !? *5 *J e ma1n rout1ne MGTEM. The 

Jtlng point; H Is the ’integration^" stl^l^TM/^ T r " 

transient in seconds, If tmax » o o « * T ? AX , 1fi the dur «t 1 on of the 

be generated. Steadi-state sol St loSs ?«!*♦ Al^ U ( c °nyerged) point will 
figure 4 are given in append 1 tout ? TI ®! 1 ? 0 po1nts In 

TOUT and H need not be the sal T h„IJ s th L pr1ntout Interval. Note that 
IBDINT. If IB 0 INT * Q, a forward dlffI™n!?I a J ?" m ? t !) od to be used Is sot by 

«fci; larasss SB 
rwms ffrravr 

nn.n, n „ tW sys t ora Ztl SiWlSia,, 

; f T t r ;"'?; ts *™ «• i" 

r.?^,r?H!2? E ;r;?:- contro ’ ,nputs ® ts-tyr^itssw 


F,4 main combustor fuel flow 


^,7 augmentor fuel flow 


a 8 


nozzle throat area 


FVOP fan variable-geometry parameter 


CVGP 


compressor variable-geometry parameter 


(high wet power). Main combustor fuel flow w u ramrma «„ , - 

i„ . , . 1 now *r 4 15 ramped In 2 sec from 0.37 

the range? of'the fa" \m 254 M 2 llZTZ. VtZVlZZT''' 

burning Is Initiated. Augmentor fuel flow w is ramnoa *. 0 
e A . * 1 r,uw w r 7 15 ramped In 3 sec from 0 to 

” ”" e ^5* *«■£ & <;??„*?; r?«r no " ,e 

steady-state run Is requested. Thus necVtu.M tiL^VS 11 ? 1 !° 0 “ hen a 
JSS to be sot equal to 0 Not* that 4 « *1^ n V??^ * 0.0-* In DIGTEM causes 

both biased and Invertedbefore leaving TMRSP* ? Th1s i! G h and . CV6P values are- 
the map Interpolation routines ° 9 ™ RSP ‘ Th1s 1s d °ne to accommodate 
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OUTPUT - TEST CASE 


The previously defined transient Is used as a test case. The lowest power 
operating point, operating point 3, Is used as the Initial condition for the 
transient. The engine Is to be accelerated from low power to full afterburning 
In 20- sec. The main routine DIGTEM to run this transient Is shown In figure 8. 
Note that backward-difference Integration Is desired (I29INT » 1); the Integra- 
tion time step Is to be 0.01 sec (H » 0.01); (ho printout Interval Is to be 
0.1 sec (TOUT - 0.1); the transient duration Is to be 20 sec (TMAX * 20.0); 
the Initial condition Is operating point 3 (NQPER * 3); Internal logic Is to 
be used to determine when a new Jacobian matrix Is needed (IHPCNV « 0); and 
there are 16 state variables (N * 16). 

The output listing from DIGTEM Is shown In figure 9. At the Initial con- 
dition (In this case, OPERATING POINT NUMBER 3) a detailed printout of the 
engine parameters Is generated at TIME » 0.0 SECONDS. This printout corre- 
sponds to the user-supplied INPUT DATA. However, the Input values may be 
slightly different from the Input data because of the effects of the correction 
coefficient scaling described earlier. Pressures, temperatures, temperature 
derivatives, mass flows, mass flow derivatives, stored mass, stored mass 
derivatives, energy derivatives, and enthalpies are listed for the various 
engine stations. Below the table, low and high spool speeds and their corre- 
sponding derivatives are listed along with main combustor and augmentor fuel 
flows, bleed flows, and variable geometries. The variables F SHIFT and CSHIFT 
are near zero since the specified values for the variable geometry positions 
result In small flow shifts out of the fan and compressor flow shift maps. 

A second printout Is generated If Implicit Integration Is selected. This 
second printout Is generated after DIGTEM converges to a balanced steady-state 
condition. If explicit Integration Is used, a startup transient will occur (If 
all the control Inputs are held constant) because of any nonzero derivatives 
that exist at the Initial operating point. When explicit Integration Is used, 
DIGTEM prints out the Input data table as described for the Implicit Integra- 
tion and then prints out a message 

FORWARD DIFFERENCE IN1EGRATI0N IS BEING USED. IF ALL THE DERIVATIVES ARE NOT 
CLOSE TO ZERO, A TRANSIENT SHOULD BE RUN TO BALANCE OUT THE ENGINE BEFORE 
CHANGING ANY OF THE CONTROLS IN TMRSP. 

Since Implicit Integration Is being used here, DIGTEM Iterates to a converged 
condition and then prints out CONVERGED STEADY STATE POINT and again gives a 
detailed printout of the engine parameters at TIME » 0.0 SECONDS. Note that 
all of the derivatives have been driven to near zero and also that the con- 
verged data are very close to the Input data. This Indicates that the Input 
data along with the correction coefficients calculated at the dry design point 
led to a nearly balanced engine at the operating point. Steady-state results 
are given In appendix E for all five DIGTEM operating points. Note that If 
explicit Integration Is used, this second printout of converged data does not 
occur. 

Next DIGTEM prints out transient results at each specified printout point. 
This Is done for both Integration schemes. Shown are TIME and the pressures at 
all engine stations In row 1; temperatures at all stations In row 2; and 
speeds and control Inputs In row 3. MATTOT Is shown In row 3 when Implicit In- 
tegration Is used. This Is the total number of Jacobian matrices calculated to 
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far\ 9 Q°lt 1n A* h JiJ| r llS 1 IJ t ;h. 1 ! *J 1s , c ? Se d ? ta are pr1 " ted out ever V 0 *1 sec 

ntPTFM At the end of the desired transient run (TMAX * 20.0 In this case) 
a K " pr nts ? u 5 a data J led ]1 §t of the engine parameters at all sta- 
^!'ni T 5 ? ccurs ^dependent of the Integration scheme selected. Figure 10 

low 0f t S e I esu1ts ' shown are Plots of high rotor speed N H 

low rotor speed Nl, combustor pressure Pa, turbine Inlet temperature T a n 

combuJ?!Ir n f 0r i t fi (>eratUre ? s funct l° ns time. As shown In figure V, main 
?SLl ?! IJ° W wa * r ? fflped f °r 2 sec fl nd then held constant. All five varl- 

S5UM « saw ‘TJFZTMriU, 

’iSrUjff t2 U ^s B f1«r5a^e? 9 except the a “ 9roent0r “" pera,ure ’ 


INTEGRATION TIME STEP STUDY 

t1me step for the test case was 0.01 sec. The 20-sec 
a^/al™* took , 13,1 k sec ° f central processing unit (CPU) time on the IBM 
ml+rw? com P u t er when the Implicit Integration scheme was used. Six Jacobian 
™ s were generated during the transient. To determine the effect of time 

warierf^f rnm n JS!**' 1on re sponse and on CPU times, the Integration time step was 
^hn^+hl 0 ?*?* 0 ? 1 ! e ?« to i; 0 sec for the same 20-sec transient. Figure 11(a) 

IteD oa 3 a f ?!? t rt ? f rmi m ? 4 Step ° n CP !J For th6 °- 001 -sec Integration time 

!!?£^?? ,3 J Se ? °f^ CP V t me was needed for the 20-sec transient. This Is 

k ?. th# I ar9 ? T5? r of steps (1 ‘ e - passes through the model). 
Increasing the- time step to 0.01 sec caused an 87-percent- reduction In CPU time 

!?r«!onf ?? C rmi A t urt ? er JH? rease to 0.1 sec caused another reduction of 69 

.Mj’Sriirtto 3?69 S e'" a11y " ’ nCrease t0 10 *« cause<1 another 

f^! 9 oS e 11( ? ) shows the corresponding number of Jacobian matrices needed 
t !’ ans1en J w1tp the various Integration time steps. As the time 
?!?? Increased, the number of Jacobian matrices needed for convergence also 
Increased from two matrices at the 0.001-sec time step to 31 at 1.0 sec. For 
6SS *£" °‘?1 sec the CPU time was primarily a function of the 
number of passes through the model. For larger time steps, however, the 

L n ? r ? d a ?? b1an !? at *l ce ? (an ? subsequent Inverses) contributed a great 
deal to the CPU time and offset much of the expected speedup. 

, l n a11 cases, stable and converged solutions were obtained. However, with 
the 1.0-sec time step, some problems occurred early In the transient with 
Inputs to maps and functions going out of range. Also some damped oscillations 
J ? [I S??®I v ® d ; K F ; r t1dle s ta Ps between 0.001 and 0.1 there was little dlffer- 

uHh fh! n Si t ?L t i ra !! S l? nt r ? sponses ’ figure 12 shows a comparison obtained 
wifi! 1,0 J ffle steps. Low rotor speed Is shown In figure 12(a). 

at T !l° sec th ®r® was a large speed difference between the two 
^m.+?J SeS ’ J p 1s occurred because the Inputs to the map and function Input 

+£?! T? °? t: , 0 n range> By 2 *° sec the H - 1 -0-sec response recovered but 
valJe of ?h° sec . and f1na11y ^owed some oscillations about the final 

?hl ? 4 ?I* S ! d * I h ?u C0 ? bu * t ? r pres sure responses are shown In figure 12(b) for 
the first 2 sec of the transient. Note the loss In accuracy for the larger 

alr^Mn! 0 A1t 5 ough no I shown In figure 12, the high rotor speeds 
and turbine Inlet temperatures exhibited the same characteristics. 
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R „ h niJl? a J! y n?Jf2 ns1 ? nt !?!V un “f 1 !!? the explicit Euler Integration scheme, 
supplied In DIGTEM. To obtain a stable solution, the Integration time step had 

;9 So.. 1 !. 8 than or ! £, ! a1 t0 0,1 m@c ' For the sai1l,e 20-sec transient, 417 sec 
of CPU time was needed. . . 


SIMULATION Of OTHER CONFIGURATIONS 

contains normalized component maps and a generalized aerothermo- 
dynamlc treatment of Its components. It also can scale the analytical model to 
match a user-specified design point. These features make It useful for slmu- 
latlng turbofan engines other than the one described In DIGTEM. Also, with 
minimal fortran reprogramming, variations from a turbofan engine such as a 
turbojet or turboshaft engine can be simulated. With major modifications to 
the coding It Is possible to model arbitrary engine configurations. 

To simulate an engine such as a turboshaft, the user need only mask 
(comment) out those areas of code that are not needed and equate variables 
where needed. The order of the state variables has been set to facilitate the 
required modifications to the Implicit Integration routine. Simulation of a 
turbosbaft engine Is described In appendix F. 

For particular engine configurations some change to the state variable 
order may be necessary. The user Is cautioned that the state variable deriva- 
tives must also be ordered as described In the section MODEL DESCRIPT-ION. For 
example, one may wish to simulate a single-spool turbojet such as the one shown 
r I. comparing this configuration with the turbofan configuration 

of figure 2, It Is clear that the fan duct, fan, and low-pressure turbine must 

be eliminated. A suitable state variable and state variable derivative order 
Is 


VS(1) « XNH 

VD0T(1) a OXNH 

VS(2) « W3 

VD0T(2) = 0W3 

VS(3) a T3 

VD0T(3) = DT3 

VS(4) a W4 

VD0T(4) a DW4 

VS( 5) « T4 

VD0T( 5) a DT4 

VS(6) = W6 

VD0T(6) » 0W6 

V$(7) a T6 

VDOT (7) = DT6 

VS(8) a W7 

VDOT ( 8) » DW7 

VS(9) = T7 

VDOT (9) » DT7 

VS(10) a WG6 

VDOT(IO) a DWG6 


Variables must be eliminated or equated as follows: 

“BUT ■ 0 03) 

FVGP a CVGP a 0 M 


W 
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a 0 


(15) 


W 2.2 



( 16 ) 


14 


p 2.2 

" p 2 

07) 

t 2.2 

- 

<1B) 

p 6 - 

p 4.1 

(19) 

t 6 " 

t 4.1 

(20) 

*6 ‘ 

w 4.1 

(21) 


In the main routine OIGTEM the number of state variables must be reduced 
to n ■ 10. The Fortran recoding to accomplish the variable changes Is done In 
the OSGNPT and the appropriate ENG1 subroutines. The state variable reordering 
Is done In the ENGBD, TMRSP, and BDPRN7 subroutines. 


CONCLUDING REMARKS 

The design and development of aircraft propulsion systems depends to a 
large extent on computer simulations. The generalized computer codes for 
developing these simulations must be flexible In being able to model many dif- 
ferent engine configurations and also must predict engine performance In both 
steady-state and transient operation. Once an engine configuration Is picked, 
the simulation must then model the specifics of that engine. Generalized 
codes, however, do not lend themselves well to modeling specific engines 
because of their generality. Also, the simulations must perform the engine 
calculations In a reasonable amount of computer time. 

Until now, the generalized computer codes available performed some but 
not all of the above functions, 0IG1EM, the digital turbofan engine model 
computer code presented In this report, has been shown to provide all of these 
functions. Besides being able to model many different configurations, DIGTEM 
provides both steady-state and transient capability, and scales Itself to match 
engine operating-point data and thus tailors Itself to model specific engines. 

DIG1EM provides all of this capability at the expense of requiring much 
more user Interaction than the other generalized codes. However, It Is written. 
In such a manner that even someone unfamiliar with gas turbine engine simula- 
tions can modify and use the simulation. To do so does require the user to 
have knowledge of Fortran. OIGTEM contains both implicit and explicit numer- 
ical Integration schemes. It Is segmented on a component basis (each component 
and mixing volume Is In Its own subroutine). Thus It can be used to do numer- 
ical Integration studies using Integration methods other than those supplied 
with the computer code. Also, because of the segmentation, parallel processing 
methods can be studied. Open-loop control Implementation Is described In 
DIGTEM. Closed-loop controls can be Implemented by adding control equations 
and Integrating the controls and state variables simultaneously or by using 
the subroutines In appendix C of Sellers, which were derived to be compatible 
with the modified Euler solution method. In addition to being a useful tool 
for simulation research and development, DIGTEM provides the flexibility to 
study a variety of engine dynamics and controls problems. 
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APPENDIX A 
SYMBOLS 

c C2?5~5 ect1onal area, cm* (In*) 
altitude, m (ft) ; 

honle flow coefficient 
nozzle velocity coefficient 
correction coefficient 

variable-geometry parameter, deg 

thrust, N (ibf) 

fan variable-geometry parameter, deg 
functional relationship, 1 =>1,30 
fuel-to-alr ratio 

COnStMt - ,0 ° “ K,/N soc 2 (386.3 lbm W1W see?) 

5?!* 1 ;? vala ® ? f fu el, J/kg (Btu/lbm) 

JZlVV c e D thal PV» J /kg (Btu/lbm) 

enthalpy change, J/kg (Btu/lbm) 

turbine enthalpy drop, J/kg *1/2 rpm (Btu/lbm °pl/2 - n(B \ 

SS.SMJB!“j Kt/iSH ”, , 

a “?W2 0 s r e c P 5^S 4 ^° S 2 R 5 0eff,c1ent > » 2 «c2/k,$ cm* K n/BtU> 

coeff1c 1ent, N 2 ,«2/k # 2 cm* K 

f^h l0n °l * l9h ~Pressure-turb1ne bleed doing work 
° f lo *“P ressu re-turb1tie bleed do1ng 9 work 
?»*“•") ,0 “ coefF ’ c,ent ' 1,2 sec2 /kg2 9 C m* k (Ibf 2 sec 2/ 

U°“ s fhTcm r (iJ U J b ’ ne d ’'' Cf,ar9e l>ressure 1055 coefficient 

Mach number 

rotational speed, rpm 

total pressure, N/cm* (psla) 

pressure ratio 

static pressure, N/cm* (psla) 
torque, N cm (In- Ibf) 

?«f a ? 0 ? Stant \ N cm/kg K < 1n Ibf/lbm »R) 
total temperature, K (°R) 1 

temperature ratio 
temperature rise parameter 
time, sec 

Internal energy. J/kg (Btu/lbm) 
volume, cm 3 (in 3 ) 
velocity, cm/sec (In/sec) 
stored mass, kg (lbm) 

mass flow rate, kg/sec (lbm/sec) 
corrected mass flow rate, kg/sec (lbm/sec) 
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Wp turbine flow parameter, kg K cmZ/N rpm sec (ibm °R In /Ibf rpm sec) 

X,Y map Inputs 

x variable 

Z map output 

0 Interpolation constant 

4 ratio of total pressure to sea-level pressure 

Y ratio of specific heats 

n efficiency 

o ratio of total temperature to standard-day temperature 


Subscripts (note that subscripts may be combined, e.g., w P 

F *4' 


A 

air 

. 

AB 

augmentor 


a 

actual value 


am 

ambient 


B 

main combustor 

r 

BL 

bleed 


BLHT 

high-pressure-turb-lne cooling bleed 


BLLT 

low-pressure-turbine cooling bleed 


BLOV 

overboard bleed 


C 

compressor 


cr 

critical flow 


0 

design Input 

, 

! 

E 

exit nozzle plane 

. 

es 

expelled nozzle shock 

1 >> 

F 

fuel 


fan 

fan 

: 

* 

i 

H 

high-pressure spool 


HT 

high-pressure turbine 


I 

Inlet 


1 

Initial condition 

\ 

10 

fan hub region 


Id 

Ideal 


In 

Into volume 


J 

station, j - 0, 2, 2.1, 2.2, 3.4, 4.15, 6, 7, 6, 9, 13, 16 

; 

y 

entrance to volume at station j. J » a, 7, ia 

? i 

L 

low-pressure spool 

i 

LT 

low-pressure turbine 

i 

load 

load 

i 

M 

map 


n 

net 

( 

new 

new 


00 

fan tip region 


old 

old 


out 

out of volume 


X 

upstream side of shock 


y 

downstream side of shock 




Superscripts: 


( )* sonic flow condition 

( * ) derivative 

-1 Inverse matrix 

Computer variables: 

AE nozzle exit area, cm 2 (In 2 ) 

ALT altitude, m (ft) 

AQL6 augmentor reactance, kg cm 2 /M sec 2 (Ibm 1n 2 /1bf sec 2 ) 

AQL13 duct reactance, kg cm 2 /N sec 2 (lbm 1n 2 /lbf sec 2 ) 

A8 nozzle throat area, cm 2 (In 2 ) 

BSCV6P bias on CVGP, deg 

BSFVGP bias on FVGP, deg 

CON nozzle flow coefficient 

CSHIFT • -change In compressor flow due to variable geometry 

CVGP compressor variable-geometry parameter, deg 

CVN nozzle velocity coefficient 

PELT change In time, sec 

OELTAV vector change In guess variable 

0H4 enthalpy change across high-pressure turbine, j/kg (Btu/lbm) 

0H41 enthalpy change across low-pressure turbine, J/kg (Btu/lbm) 

0T3 temperature derivative In compressor mixing volume, deg/sec 

0T4 temperature derivative In combustor mixing volume, deg/sec 

0T41 temperature derivative In high-pressure-turbine mixing volume, 

0T6 temperature derivative In low-pressure-turbine volume, deg/sec 

0T7 temperature derivative In augmentor mixing volume, deg/sec 

0T13 temperature derivative In fan mixing volume. deg/seC 

0WA13 fluid momentum derivative In duct, kg/sec 2 (Ibm/sec 2 ) 

0W3 flow derivative In compressor mixing volume, kg/sec (Ibm/sec) 

0W4 flow derivative In combustor mixing volume, kg/sec (Ibm/sec) 

0W41 flow derivative In high-pressure-turbine mixing volume, kg/sec 

(Ibm/sec) 

0W6 flow derivative In low-pressure-turbine mixing volume, kg/sec 
(Ibm/sec) 

0W7 flow derivative In augmentor mixing volume, kg/sec (Ibm/sec) 
0W13 flow derivative- In fan mixing volume, kg/sec (Ibm/sec) 

DWG6 fluid momentum derivative In augmentor, kg/sec 2 (Ibm/sec 2 ) 

OXNH high-rotor-speed derivative, rpm/sec 

DXNl low- rotor-speed derivative, rpm/sec 

E error vector 

EMAT Jacobian matrix 

ERRBSE vector of past errors 

ETAAB augmentor efficiency 

ETAB combustor efficiency 

ETAHC compressor efficiency 

ETAIF fan hub efficiency 

ETAOF fan tip efficiency 

FG gross thrust, N (ibf) 

FN net thrust, N (Ibf) 

FRAC external control for matrix convergence 

FVGP fan variable-geometry parameter, deg 
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FSHIFT 

H 

IBDINT 

IHPCNV 

IPRINT 

ISS 

JSS 

MAPNO 

MATRIX 

MPAS 

N 

NCOM 

NCV 

NFCT 

NOBUG 

NOPER 

NPT 

PO 

P2 

PI 3 

P22 

P3 

P4 

P41 

P5 

P6 

P7 

PCNCHG 

POINT 

REF 

TAM 

T2 

T22 

T3 

T4 

T41 

T6 

T7 

T13 

TMAX 

TOL1 

TOL2 

TOLPCG 

TOLSS 

TOUT 

V3 

VI 3 

V4 

V41 

V6 

V7 

VDELT A 

VOOT 

VDOTSV 

VDOTT 


change In fan flow due to variable geometry 

time step, sec 

Integration select switch 

matrix calculation select switch 

print select switch 

initial-condition switch 

steady-state switch 

Indicator for component map 

switch for generating a new Jacobian matrix 

maximum allowable Iteration passes 

system order 

map Interpolation switch 


number of curves on a map 

number of common functions 

debug printout select switch 

operating-point select switch 

number of points on a curve 

Inlet pressure, N/cm 2 (psla) 

fan Inlet pressure, N/cm 2 (psla) 

duct pressure, N/cm 2 (psla) 

compressor Inlet pressure, N/cm 2 (psla) 

combustor pressure, N/cm 2 (psla) 

high-pressure-turbine Inlet pressure, N/cm 2 (psla) 

low-pressure-turbine Inlet pressure, N/cm 2 (psla) 

Intermediate pressure, N/cm 2 (psla) 

augmentor Inlet pressure, N/cm 2 (psla) 

nozzle Inlet pressure, N/cm 2 (psla) 

Iteration convergence rate 
operating point 

desired value of summation of errors 

ambient temperature, K (°R) 

fan Inlet total temperature, * (®R) 

compressor Inlet temperature, K (°R) 

combustor Inlet temperature, K ( °R) 

high-pressure-turbine Inlet temperature, K (*R) 

low-pressure-turbine Inlet temperature, K (*R) 

augmentor Inlet temperature, K (*R) 

nozzle Inlet temperature, K (°R) 

duct temperature, K (°R) 

transient duration, sec 

lower limit for partial derivative 

upper limit for partial derivative 

convergence rate at which a new Jacobian matrix Is generated 

error tolerance 

output time step, sec 

compressor volume, cm 3 (In 3 ) 

duct volume, cm 3 (In 3 ) 

combustor volume, cm 3 (In 3 ) 

high-pressure-turbine volume, cm 3 (In 3 ) 

low-pressure-turbine volume, cm 3 (In 3 ) 

augmentor volume, cm 3 (In 3 ) 

Initial perturbation In state variables for matrix generation 
vector of state variable derivatives at current time 
vector of state variable derivatives at previous time 
vector of average state variable derivatives 
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WA13 duct flu?d momentum, 6g/sec 2 (ibm/sec 2 ) 

WA2 mass flow rate at station 2, kg/sec (Ibm/sec) 

WA22 mass flow rate at compressor Inlet, kg/sec (Ibm/sec) 

WA3 mass flow rate at combustor Inlet, kg/sec (Ibm/sec) 

WG4 mass flow rate at high-pressure-turbine Inlet, kg/sec (Ibm/sec) 

WG41 mass flow rate at low-pressure-turbine Inlet, kg/sec (Ibm/sec) 

WG7 mass flow rate at nozzle, kg/sec (ibm/sec) 

W13 duct volume stored mass, kg (Ibm) 

W3 compressor volume stored mass, kg (Ibm) 

W4 combustor volume stored mass, kg (Ibm) 

W41 high-pressure-turbine volume stored mass, kg (Ibm) 

W6 low-pressure-turbine volume stored mass, kg (Ibm) 

W7 augmentor volume stored mass, kg (Ibm) 

WBLHT high-pressure-turbine cooling bleed flow, kg/sec (Ibm/sec) 

W6LLT low-pressure-turbine cooling bleed flow, kg/sec (Ibm/sec) 

WBLOV overboard bleed, kg/sec (ibm/sec) 

WF4 main combustor fuel flow, kg/sec (Ibm/sec) 

WF7 augmentor fuel flow, kg/sec (Ibm/sec) 

WG6 augmentor fluid momentum, kg/sec 2 (Ibm/sec 2 ) 

XIH high rotor moment of Inertia, N cm/sec 2 (lbf In/sec 2 ) 

XIL low rotor moment of Inertia, N cm/sec 2 (lbf In/sec 2 ) 

XMN Mach number 

XNH low rotor speed, rpm 

XNL high rotor speed, rpm 

XXX summation of squares of changes In errors to maximum error 
YYY state change vector 


♦* 
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APPENDIX B 
ANALYTICAL MODEL 


The mathematical model describing the two-spool, two-stream turbofan engine 
In DIGTEM Is described In detail In reference 8. Overall performance maps are 
used to provide the steady-state representations of the engine's rotating com- 
ponents. Fluid momentum In the bypass duct and the augmentor, mass and energy 
storage within control volumes, and rotor Inertias are also Included to provide 
transient capability. For completeness, the mathematical model Is presented 
below. 


Steady-State Model 


Flight cond i tion and Inlet . - The following conditions define the flight 
conditions and Inlet model: 

Pq ° D( a ) 


To ° f2( a ) ♦ T am 

m » i .0 if Mq < i.o 

- 1.0 - 0.075 (M 0 - I.O) 1 * 35 If Mq. > 1.0 


Wo 


1.0 + 


(yt - m 


(Bl) 
(B2) 

CB3) 
(B4) 

(B5) 
(B6) 

where functions fi and f 2 are curve fits to atmospheric data from 
reference 10. 

Gas properties . - Curve fits of data found In reference 11 are used to 
compute variable thermodynamic gas properties. JP-4 Is assumed to be the fuel. 
do£~each- .control volume- the following equations are used: 



« f 3 (T,f/a) 

(67) 

f4(f/a>«R A 

(B8) 

• c D -| 
P J 

(B9) 

c. 

Y * 

1 c 

(BIO) 

V 

- fs(T,f/a) 

(Bll) 
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performance Is represented by a set of ovorall performance maps. 

reS^n+ m JEn JUfnrm d for Jt!! t1p a a !!? hub se <: t1ons ‘ The mapsare assumed to 
♦?SIJ 8en i.! a ! 5! rf 2 r ?"* e w1th var1ab1e 0 eom0tr y at nominal, scheduled posl- 
tlons. Map-generated, fan-corrected airflow Is adjusted to account for off- 

schedule geometry effects. The following equations describe the fan model * 


OfflGINM- W* *} 
OF POOR QUAUTt 


(wJ 


fan,M 

P 2.1 " P 2.2 ’ P 2 f 

■ 

(w c ) i 7 1 + f 

a . kn-H li 


(hi UlA 

\Y.<7 

, (hi !L_\ 

’W*l 7 

fem, N L /eJ /2 l 


e 


1/2 


n fan,0P “ f 9 




'1.,.,. ■ ft) 




fan 


( 812 ) 
(813) 

(614) 
(815) 

- 1.0 (B16) 

(817) 

(818) 

- 1.0 (B19) 

! (B20) 

(821) 

«h 4 i S R L ?^°r , >: °m al ] P® rfo ™ante ™P* ar ® "*** for the compressor with a 
shift In the corrected airflow based on off-schedule values of variable- 

geometry position. The following equations describe the compressor model: 
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] t 2.; 


fiir«h?if^ 7 ‘Ai 4 H! low ^ throu ? h *5®. ? l6 ! d P asSa 9®s Is assumed to be choked. Botl 
tjjblne cooling and overboard bleeds are modeled. The equations are as 


fe) 6L - p 3 (K) (^rh) 


(Y 3 +1)/2(Y 3 -1) 


W BLHT “ A BLHT Ia 


W BLLT * A BLLT (a 


“blov “ a blov Va 


BilMnis. - Overall performance of the hlghand low-pressure turbines Is 
represented by bivariate maps. Cooling bleed for each turbine Is assumed to 

toi-blno discharge although a *rttJTo» 
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(037) 

(B38) 

(039) 


Combustors and ducts . - Total pressure losses are Included In the models 
of the main combustor, bypass duct, mixer entrance, and augmentor. Heat addi- 
tion associated with the burning of fuel In the main combustor and augmentor 
Is assumed to take place In volumes Va and V 7 , respectively. The following 
equations describe the combustor and duct models: 

1/2 


• a 

w 3 


p 3< p 3 - p a> 
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K B T 3 


V.?bV ( ] - V T 4 
Oh B « HVfn B 

n B - f 10 t(f/a) 4 ] 
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P ) 
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p « - 


( B40) 
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T AB " B AB T 6 + (1 " »AB> T 7 

(B47) 

Ah Afi " HVFi, AB 

(BAB) 

n AB B 

(B49) 

Wr 7 t Wp j 

(f/ 4) . SJ-—IA 

W 6- W F,4 

(B50) 

f . p » ' K o“i3 T ia 
6 P 13 

(B51) 

T 6 ‘ T 13 

(B52) 


nozzle. - A convergent-divergent nozzle configuration Is assumed. 
The following equations define the basic nozzle model and are based on material 
from reference 12. Simplifications to the basic model, Intended to reduce 
computation time, are noted: 
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If Pq/P? > (Pq/P-*) » the flow Is subsonic In the nozzle and 

cr 
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Otherwise a shock may exist In the divergent portion of the nozzle. To compute 
the required parameters under these conditions , shock tables such as those In 
reference 12 must be used. . 
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If p o/ p 7 • ( r o/ p ?)es* s.hock will be In the nozzle exit plane. Then 



(B68) 
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If P 0 /P ? < (lytyes* the shock 1s e * ternal to the nozzle. Then 
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> *V P 7 > ^ P 0^ P 7^ * the shock 1n the divergent section and 
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To solve these equations, M x can be varied until equations (882) and (877) 
produce the same values for Pg. Then 
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mK /r t1ons ( B53 > to < B83 ) are 1n subroutine NOZZL. The Inputs to nozzl aro 
ambient pressure, nozzle Inlet total pressure, nozzle Inlet total t!!mperatu% 
nozzle throat area, nozzle exit area, nozzle pressure ratio, nozzle flow coef* 
f Iclent , and nozzle velocity coefficient. FUN1 Is called b^ NOZZL to 

terpolate tabular data of functions f^Aj/A,), fj| (P./P,). f 

f24( A E^8)» ^nd f 30 ( A^/Aq) . Functions f 25 <Mw) and f 9 c(M v ) are reoresenterf Kw 
quadratic functions. The Iterative loop associated with’ shSck InthedlverSent 
section Is replaced by a quadratic function of prossure ratio and is kiacoh^ku 

exlt^veloclty; 00 ° f 4r " rat ’°- rasu,t * M * 


V E ’ M E C v 


fa»h*A U2 

■ H \ 'H * V 


(884) 


The net thrust Is computed by subtracting Inlet ram drag from gross thrust: 

1/2 


F n - M q w 2 


Wo 1 


9, 


(B8S) 


Engine Dynamics 


Intercomponent volumes. - As shown In figure 2, Intercomoonent volume 
assumed at engine locations where (1) gas dynamics are considered Important or 
tn +h 3S dynam1cs a f e re< iu1red to avoid an Iterative solution of the equations 

Jftfw S fh V °] Ume V tor !! 9 ? of ® ner 9y and mass occurs. The following equations 
define the dynamic models of the Intercomponent volumes (fig. .2): 


f t 

W 13 ’ )„ <*2 - *2.2 - “,3 > dt * “,3.1 

T 13 ' £ |[<*Z - *2.2) ("13 - h 13> /*».13 

* T 13 ( *2 ' *2.2 - *13><Tl3 - >)] j W, 3 j dt * T^, 

p . R A W 13 T 13 
13 »,» 
vat 

W 3 " )o ( “ 2 - 2 ' W BLH1 " “but “ w BL0V “ “ 3 ) dt f W 3 j 


( 886 ) 

(887) 

( 888 ) 
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io ’ h 3 ,/Cv - 3 * T 3^2.2 " W BLHT ' “BUT * “b10V * “3* 

«(T 3 -n]/M 3 |<it* T 3i , 

(B90) 


P 3’ V 3 

( B91 ) 

“4- j 
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T 4 ’ { (j[ 4 3 h B ‘ “f,4 th B - h 4 (i 3 * 4 F,4>] / C V, 4 

* T 4< 4 3 * “F,4 - “4 >(1, 4 - 11 j/ “ 4 ) dt * T 4.1 
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p, „ _0_a 
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W 7 ■ J* < 4 6 * 4 F,7 - 4 7> dt * W 7.1 

jj (|[“6 h AB * W F,7 lh A6 * h 7 (4 6 * 4 F,7>] j c »,7 
* T 7< 4 6 * if >7 - “ 7 )(T 7 - 1)| / W,^dt 
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(8103) 


Fluid momentum . - The effects of fluid momentum are considered In the 
bypass duct and augmentor duct models: 


w 13 * 9 c 
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6 - V dt + *13,1 
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Rotor Inertias ♦ - Rotor speeds are computed from dynamic forms of the 
angular momentum equations: 


N L “(ir) I L ( Q j[ Ah LT (w 4.1 + K BLWLT W BLLT^ " * W 2 " W 2.2 )(h 13 “ h 2 ) 

" *2.2 (h 2.2 " h 2 } ] j H tj dt + N L,1 

17 | [ Ah HT< W 4 ♦” K BIWHT W BLHT> 

- w 2 2 (h^ - h 2 2 )j /n h J dt ♦ N h>1 
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Correction Coefficients for "Trimming" Model 


In DIQTEM, design-point data throughout the engine are specified as Input. 
If the turbofan engine model In DIQTEM was exact, the specified Input data 
would lead to a perfectly balanced engine condition. However, Incompati- 
bilities between the DIQTEM Input data and the model will result In nonzero 
derivatives or mismatches between predicted and specified outputs of component 
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Jo compensate for these differences, a "self-trimming" feature has been 
hai n *° f * or [ e ^^ on coefficients are calculated In subroutine DSGNPT 

include^? thB p ng1ne T at th ® dry design point, for example, the Input data 

*2i D : P2 i°i I 2 »°\ Ta 5!: arul a D» where subscript D Indicates 

design-point Input data. In OSGNPT, subroutine FLCQND Is called by using an 

Mn* Slid as Innnfc Tho roc i » 1 4* 4 nn AuiK.ii ^ 


ng, and T am as Inputs. The resulting output variables are Po a , T? =,, 
value 3 Ideally* Wher ® subscr1pt a stands for the actual calculated’ 


P 2,a “ P 2,D 


(BIOS) 


2, a 2,0 (B109) 

P 0,a * P 0,0 < B1 1 0) 

However, If they are not equal, the equations that use these values will be 
by correction coefficients. OSGNPT Is called only once (at the design 
Sa!? 1 ^l? 1 * * bese coefficients. The correction coefficients are then 

t 5k mode ]: They are used at both the design point and the off-design 
points. The scaling coefficients for the Inlet conditions are * 

CC, « ^*1 (Bill) 


(B112) 


(8113) 


These coefficients are used to scale the Inlet model. Equation (Bl) becomes 


P Q = f^a) x CC, 


( Bl 1 4) 


equation (B4) becomes 
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J 2 ’ T o ’ •» * “4 8 


and equation (B5) becomes 
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The other correction coefficients and their corresponding "trimmed" equation 
are presented below. For the fan M 
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and equation (B14) becomes 
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and equation (B13) becomes 
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and equation (B20) becomes 
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For the compressor the same scaling procedure Is used. That Is 
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and equation ( B23) becomes 
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and equation (827) becomes 


For the turbines 
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Equations (B33) and (B37) become, respectively 
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The next set of correction coefficients zeros out the state variable 
derivatives associated with energy balances In the Intercomponent volumes: 


CC 10 a h 4.0 (*3.0 4 “r t 4.pV *3.D h B.I 


W F,4,0 Ah B,D 
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Note that 
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In steady state. Also 
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and equation (B96) becomes 
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and equation (B99) becomes 
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" M 14 Ah LT (w 4.1 + W BLIT K BUM.T*] j c v,b 
* T 6 <*4.1 * “BUT * 4 13 - 4 6> <*6 " ’>j /»(,) « * T 6j , 
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In steady state 


w 4.1 * W BLLT * w 13 ” w 6 " 0 


(B141 ) 


The next two correction coefficients are used to zero the speed derivatives at 
the design point. For the high rotor speed 
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w 2.2,0 (h 3.0 “ h 2.2.0 ) 


15 • i 

Ah HTD * W 4,D + K BLWHTQ w BLHTD* 

and equation (B107) becomes 
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For the low rotor speed 

CC 16 » (W -^. P ~ W 2.2,0 )(h 13 t D ~ h 2.D^ * W 2.2.0 (h 2.2.0 ~ h 2.D ) (B144) 

Ah LT0 <Vl,D * K 3LWLT0 w 6LLTD ) 
and equation (B106) becomes 

N l * («") T[ jj j [ Sh LT (w 4.1 * K BLWLT w BILT ) X - C il6 

- (W 2 - W 2 2 ) (hj 3 - h 2 ) - w 2 j, (hj 2 - hj)j ^N L | at ♦ N Li , (BUS) 

The last three correction coefficients compensate for the Imbalances In the 
augmentor.. and nozzle modeU 
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and equation (B53) becomes 
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In the augmentor 
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Finally for the thrust 
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F N,a ~ A E < P E " V 


and equation (B54) becomes 
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APPENDIX C 


INTEGRATION AND ITERATION SCHEMES 


▼ It 17 


* m v Mi my Ultnmquo 


to calculato'stoaSy-stat^oporatfno'point^ 1t fho t c V ? m ? t ?? cl that OISTEM uses 
n at i? 9 p f 1nt "W'ro, tSo S olu,?oS of a ' s of a s toady-stato 

responding to various engine matchlna coh'ifr>/Hn+? ?f nonlinear equations, cor~ 
airflows, compressor and turbine work funf-HAnt^ s !J cb as rotational speeds, 

To satisfy those constraints, the“ l,t be " m, ° f ’°“ functions, 
parameters that can bo varied fsach as ramnan.f ublo an equal number of engine 
and flow functions) for the tnrhnfan compressor and turbine pressure ratios 

onglneparamoters; independent ‘variables )°a!!d K 0I6 I EH thm «' 

variables). 0IGTE« seafches fo? thfvl L oJ Ib error varlables (dependent 
results in the engine error vM, ^ that 

ables by E, the P match1nrequ^ and the de P e " de "t varl- 

Ei (VSj) o 0 1„ 1,2,..., N; J a 1,2,..., N (Cl) 

The procedure used to satlsfv th* «*,♦ n « 

able Newton* Raphson method, where changes 1n 11 ? ea IeI q i at1on 2 J s the roul t1var1- 
changes In vs b, a flrstW? ffemerencTeS *° * re ' ate<l *° 


AE « EMAT x AVS 
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reference^ondltlon (oplratln^^ in M VS and E from some 

partial derivatives of E w1th P °espect to^vs ** 30 NXN Jacob1an matrix of 


EMAT 
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1 
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(C3) 


that equatlon^cl) 1? app?2x1i!iSted [[y fferences about an grating point such 


EMAT 


AE. 


U “ AVSj 1 “ 1,2 Hi j a 1,2 N 


(C4) 


operating poInM^Improve^by 15 obta1ned » the steady-state balance at the 


vs new - VS old - EMAT-1 x E o1d 


(C5) 
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If the system af equations were linear, the process would lead to convergence 
In one Iteration. In practice, nonllnearltles In the system prevent Immediate 
convergence. In this case the new VS and E are taken to be the reference 
values and a new matrix Is generated. If the system Is not too nonlinear and 
the Initial guesses for VS are reasonably accurate, convergence Is achieved 
In relatively few Iterations. 


Dynamic Equations 

Once an Initial steady-state solution has been obtained, a time-varying 
solution may be generated. This requires the solution of a set of differential 
equations that model the system. In this section the procedure used to solve 
the set of differential equations In DIGTEM Is discussed. 

Consider the differential equation 

x « f(x,t) (C6) 

lo obtain the numerical solution on a digital computer, the differential equa- 
tion must be approximated by a difference equation. One common method Is to 
use Euler's method where equation (C6) Is approximated by 

*n+l ° ^n + f (Xn» *n) *T (C7) 

Equation (C7) allows for explicit calculation of X n+ i as a function of the 
previous values of X n and t n . This Euler Integration method Is the 
forward-difference Integration scheme Included In DIGTEM. When an explicit 
method Is used for Integrating a system of equations, the Integration time step 
Is restricted by the highest frequency In the system (as derived In ref. 6). 
However, dynamic engine simulations contain both high and low frequencies. The 
high frequencies result from the lumped-volume representation of component 
dynamics, which Includes storage of mass and energy. Low frequencies result 
from rotor dynamics. In DIGTEM the range of frequencies for the test case Is 
0.4 to 330 Hz. Frequently the simulation user Is Interested In low-frequency 
effects such as rotor transients and Is not concerned about high-frequency 
effects. These transients typically are 5 to 10 sec In duration. However, the 
user must still use a small Integration time step to Insure numerical stabil- 
ity. Although It gives very accurate results, this requirement can cause large 
amounts of computer time to be used. In OIGTEM the largest time step that can . 
be used with the Euler Integration method Is approximately 0.1 msec for the 
test case. Thus the Z0~see transient used In the DIGTEM test case consumed •• 

417 sec of CPU time on the IBM 370/3033 computer. 

Another method for approximating equation (C6) Is the Improved Euler: 

X h+ i » X n ♦ f(X n+1 , t n+1 ) AT (C8) 

In general equation (C8) cannot be solved explicitly for X n *i because of the 
dependence of f on X n+ i . Thus some form of Iteration must be used at each 
time step. For this Implicit formulation there Is no restriction on step size 
(ref 6) to guarantee numerical stability. However, some loss In dynamic 
accuracy can occur If the step size Is too large. 
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whore the first guess at X^i is givon by 

^ntl - X n ♦ f (X n , t n ) at (C1Q) 

rla a r^sot 1 ^ f 'od ! * 0 ?h1 2 V ?» y tho°? n tc 0 rlt 1 on ^tRSJ u^WSweT'” "" 


Implementation In DIGTEM 

point 1s 0 i 6 d©?lgn a pSlntf correct1orcoofn?C? ; ' 1n9 po1 ? tt If ' the grating 
a balanced engine condition The C ! ents are c ^culated to try to force. 

and the correction coefficients r I f re calculated by using the Inout 
vector Is the derivative vector; thus 3 steady* state condition the error 
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rn a 4 J a1a Joed condition; If not/oiGTEM Iterates^ J ser >» *£• operating point 
earlier. °" techn,q “" i ’ 5 *25* 

a matr1x P of partial derfvatlves^^chSn 35 m ? tr1x EMAT »u*t be formed. EMAT Is 
changes In guess 5arl!bles ^ OIBTffl 8! I" #rror r ar1ables w”h mpe?? to 
variables (VS) and EMAT Is formed by "in.* difyssj** are we sta? ' c 


EMAT(I,J)«(E(J) - ERRBSE( J) )/QELTAV( I) 


• - (Cl?) 

and VSCt0r 


v W-r«»T 


W . .S...E.. > vs 


old 


(C13) 


UPdat !?\ t ?Ja e n s f^ e t r^’be'™ er ^'’ r V ' S •« within 

controls or Input conditions 1s^Ur?f* P I? ced l!r e ? as fol1 °ws: one of the 
tlon. for transient operation thl !" t"* * tead ^ sta te balanced condl* 
Improved Euler approximation ° r vector 1s redefined by using the 


VDOTT X OELT - (Vs 


new 


- vs ,J - 

old' 


(C14) 
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m r g ' n r"* " »*r r* ■ * 




where 


VOOTT 


eOT t VQQTSy 

2.0 


(CIS) 


and all errors are scaled by the last converged corresponding state variable. 

Note that In steady state VDQT equals VDOTSV and VS equals VS 

However, once an 1nput_or_control Is changed, one or mor^of the errors^are 
forced to be nonzero (VOOT changes). This forces a change In to satlsfv 

eJSatlSn (C13). Updat1n9 1s accom PHshed by using the al?eady generated EHAT and 

As ^new starts t0 move awa * the Initial operating point, the original 

fhP°?hfnn? a Jl 1X 1s used “ nt11 1t i no longer provides a good approximation of 
the change In errors with respect to the change In states The deri «ir»h fn 

calculate a new matrix Is defined by the user by setting TOLPCG D-IGTEM cal 
ciliated ? f Convergence rate ’ PCNCHS > ^Ing a transit" 9 A new matMx Js call 


PCNCHG < TOLPCG 


(C16) 


A new matrix Is also calculated If the maximum number of allowable oasses MPas 

lhe^ C c^uJ UMn9 an U ! rat1 ? n ‘ H.ra again MPAS Is set by tXe uJef loth 
are used t° try to minimize the number of Jacobian matrices 

lists the 9 imni/?I e 4 S + S s these calculations are time consuming. Table II 
lists the Implicit Integration parameter settings In BDINT6 These setting 

"Oden and data supplied w1th S DIGf E Sbl!t cap be chipped bj 
configuration^ ^ 0CCUr “ Uh d fferent ,nput data °r a different engine * 


nau ia ua icu lanon 


scheme' 1 converge! eVera1 fMt “ r " 8tlINTG that he,p the tntegratlon 


,^- r -t-«rbat1on calculation. - since finite differences are used to aenerat* 

tant J3C if thpl a ^ r X f th ? SlZ6S ° f the P ert urbat1ons of the states are Import 
tant. If they are too large, errors will be Introduced by the system non- 

linearities. If they are too small, the partial derivatives will be In error 
because of numerical problems (without double-precision arithmetic). 


*J hus ? mechanism has been Included In BOINTG to optimize the sires 

ra n iScKTs once * ^ 


XXX 


l Ji~l 

1 1=1 


E(I) - ERRBSE (I) 


] 


(C17) 
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my • *# ****■ m y * * ; 






for each state variable and then checking If 

T0L1..< XXX < T0L2 


(CIS) 


^and, the matrix Is considered "good * 
ilatlon. tqli - 0.001 and TQL2 • 0,01 work well for the operating 


If all XXX 's .... 

For this simulation, T0L1 
points. 

a P0lnl Hh2"xKK- g e wlVl t, 5 ot i> g^ ri»"*gg"? r ? T • ,0r * he <"'‘< a ’ parturitions ,t 

nnTWTr U W1 i n ?\: a11 in the tolerance band described above. Thus 

Sn^^aM 1 ^ t0 try t0 f ° rCe the XXX ' S wUh1 " the bar ^ 


YYY o 


EEL 

XXX 


(C19) 


for each state variable, 
band: 


REF Is defined as being the center of the tolerance 


REF 


i m t TOl 

2.0 


(C20) 


banS, t the S setVm?s l ’?i sto? e d a1 'J» ted .J! ch k tha J the xxx ' s fa » “’thin «• 

perturbations Sffi?* . bUI “ .7 al ? ? $ d I^r^nTuaE" th! 

z'tVstnVtiz out i^ s- t;i.?s2rJi!r‘ w j , 1 , ,o ysj‘2; the 

required for S uS?e5«nt J!U?r* SeSe?S?lS! ly ^ of passes 


Error Messages 


sets ilthlU e fElMan?e Par uhM iM!';*"'' a f' tuat,on "»» arise where XX* never 
message: tolerance * when th1s happens, the program prints out an error 


CHECK INPUT - BAD PARTIAL,. DERIVATIVE 

tsnxs r^sasi-.M sns& a: 


ITERATION FAILURE 15 51 20 


?r b a e nd 

the point at which the convergence failed (ITRAN). P 5 1 '» and 


In this situation, a debug output Is printed. This Is the same debug 
output a for the bad partial derivative, and It Indicates 

I counter up to system order 

VS current guess variable 

VCONV past converged guess variable 
VOOT current state derivative 

VD0T1 averaged state variable derivative between current time point and 
past time point 
E current errors 

After the printout the simulation continues. Note that with the Implicit 
method a convergence failure can. occur even If the errors are very close to 
the tolerance band. Since the simulation may recover after the failure, the 
simulation Is allowed to continue and the user may make a Judgment as to the 
validity of the data after a convergence failure. The occurrence of many con- 
vergence failures In a transient, however, usually Indicates a need for the 
user to Increase tolerance or to check the Input and coding. 

The debug output, as described, Is generated by BDINTG by setting 
NOBUG » 1 (table II). The user may want to reprogram the logic to obtain 
debug output at other times whan-dlfflcult convergence problems are encoun- 
tered. 

Other ert&Mnessages are Issued In DIGTEM. These are 
MAP. NO. INPUTS OUT OF RANGE 

XIN = YIN = MATRIX . 


and 


FUNCTION NO. INPUT OUT OF RANGE 

XIN * MATRIX « 

These are output from subroutines MOOR and FOOR, respectively. The MAP NO In 
the first error message corresponds to the MAPNO described earlier for the 
component maps. The function out-of-range problem Is a little more difficult 
to debug since the single-valued Interpolation routine FUN1 Is used In sub- 
routines FLCOND, TRAT , and NOZll. In either out-of-range case the Inputs are 
printed and the user must locate the map or function In question to debug the 
problem. Depending on the engine being simulated, maps or functions, or both, 
may have to be extended. MATRIX Is printed out to Indicate If a perturbation 
Is being performed to generate a Jacobian matrix since this may cause a map or. 
function to go out of range. 


Convergence Aids 

BDINTG has some built-in parameters to help the user If convergence 
problems occur. These are listed In table II. FRAC can be used to force a 
larger or smaller Iteration time step. T0L1 and TOL2 can be shifted depending 
on the linearity of the system being simulated. TOLSS can be Increased If 
convergence Is difficult. MPAS can be Increased and finally TOLPCG can be 
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<*)• 


2Jto? 5 S SJ« k ISS s ?t Internally to define a steady-state or transient run 
MATRIX can be controlled externally by using IHPCNV (from table I) ipI>int’is 

Mtar Mts nMii » voelta - >* 



APPENDIX D 


p' “ r,v w. 


FLOW CHARTS 


POOH ij U/iilT/ 


This appendix contains flow chajcls-fon-the-maln p/jxgram DIGTEM and all of 
Its subroutines. 


Writs DIGTEM heading 


Specify time step, printout, 
operating point, and transient data 


Read in maps and steady- k-^C ftAT A 
state design -point data L^Jinuatai 


JS1 Design No 
\ point? / 


Calculate correction 
coefficients 


JDSGNPTl 


Calculate engine 
parameters 


Set up guess vector and 
derivative vector 


Backward- 

difference 

integration? 




Do backward- 
difference integration 

















OKWHi:U. r.'C 

OF PQQK QUALltV 


/Print out Indication of 
bad partial derivative. 
Print out englno 
varlablos / 


Converged? 


Updato print counter 


Cannot get a good 
partial derivative - 
no good Jacobian 


Generate a 
Jacobian matrix- 
if possible 


[Good matrix 


Invert matrix 


Print desired, 
variables 


Update guesses 


Calculate open- 
loop controls 


Form new 
guesses 


Calculate engine 
variables 


Calculate errors 


Converged? 


Save converged 
states and 
derivatives 


fr^ ^ *”**** 1 ^ * L v 












Subroutine DSGNPT 



Calculate design 
enthalpies 



Calculate afterburning 
correction coefficients 



Equate design parameters 
to Input 


1 * 


Calculate fan Input 
total temperature 
and pressure 

^flcono) 

r ^ 

Calculate design 
enthalpies at 
various stations 


♦ 


Calculate temperature ratios 
from pressure ratios at 
various stations 


l 

Calculate nonafterburning 
correction coefficients 



(DIGTEM) 


* 

i 


\ 


\ 

» 
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QRtGliM* HACK; m 
OF POOR QUALITY 


Subroutine) DMINV 



Do a doublo-proculon 
matrix Invert 



Subroutine DUCT 
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ORIGINAL PACtt- 
OF POOR QIJALCiY 












ORIGINAL f‘W2E Jj^ 
OF POOB QUAUTY 




@c[ 

Compute low-pressure- 
turblno performance 

T 

(en^C^I 

Compute hlgh-prosturo- 
turblno mixing volumo 
dynamics 


7 

(0d 

Computo low-prossuro- 
turblno mixing volumo 
dynamics 

Qp\ 

♦ 

1 

Computo nozzle porlormaneo 

t 

Qfz 

Compute augmontor perfor- 
mance and augmentor 
mixing volumo dynamics 

^NG^f7 

^014^*7 

* 

Computo spool dynamics 

J 

Compute duct dynamics 


/calling 

routine 
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IV- ^■TT’'1'W 


ORIGft'AU t’*W '’ : - 

OF POOR QUAU (TV 


Subroutine EflCl 



Subroutine ENG2 


Calculate compressor 
map Inputs 


Find flow from map data 


Find efficiency from 
map data 


Calculate flow due to 
CVOP shift In map data 



Subroutine ENG3 




Subroutine ENG4 
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ORIGINAL PAGE IS 
OF.. POOR..QUALITY 


Subroutine FNGS 


Form continuity equation for 
compressor mixing volume 


Calculate pressure ratio and 
average temperature 


Calculate temperature ratio 


Calculate enthalpies ITiPROCOM) 


Form energy equation 


Calculate volume 
stored mass 


Perform Integration 
calculations 














Subroutine ENG8 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


Suhroutlno FMNTG 
[DIGTEM^ 



Sot translont Indicator and print counter 

I 



Set open-loop controls 


Calculate engine response 

~T— 



ENGINE 


Increment print counter 


Printout 

desired?- 


Yes 


No 


Print out variables 


]^TPRIN^ 


Reset print counter 
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Subroutlhe TPRtNT ORIGINAL PAGE m 



Subroutine TRAT 










APPENDIX E 

DI6TEM STEADY-STATE OPERATING POINTS 

th« f?IS T c?+l? as e U 4 ed 9® n ® ra J® converged steady-state operating points Tor 
s ®* 5 ,®f Input data. To obtain these data, NOPER was set to the 

DIGTEM d 0 irnTrT?M P Jh nt ^ nUm S er 4 and T ? A * was set to 0,0 1n the ma1n program 
Th«Io M H a + In DI6UM *J e dr 5f design point- Is the first operating point read In,. .. 
These data are used to calculate correction coefficients for balanclna the 

the model^t all * +I heSe sam ® ^ correction coefficients are then used'?? «ale 
!!«' l 1 t t er 0 ?? ra *! ng Points. A similar procedure Is used for wet 
operating points, where the first wet operating point read In Is the wet desian 
;; "*• c ; r [« ct 1on coefficients that only affe?t wet opJ?at Sg points lr6 cl\ 
{*• wet desjon point and applied to the other wet op£a?lS, JolS?s 
I?I«i?* ,eC 4 l on ®°® ff Sclents are Imbedded In the DIGTEM equations and serve to 
Ityut data!* n0t t0 el1ra1nate ’ d1 ^erences between the OIGTEM model and the 

, i.[ 1g ^r e 14 shows Printouts of corrected and converged steadv-state data 

J* a the ' l:%2S? r ;i 1n " po1nts ; .[’"•« «(•) contalnslhe d?y d2s 1 g^pol St 

thl e T? des,an po,nt sre «n close to zero 

oecause of the scaling by the correction coefficients. DIGTEM then Iterates 

sa’s 

Th.sis 5iiei]?s*c!:f)?;{^, , t:?" > o:{ c :f:tS , :rtK 1 ^, w 

JalJulated^rthrSffd^? off i e ! 1gn points as well. Correction coefficients 
calculated at the off-design points would be slightly different than those 

calculated at the design point because of slight Inconsistencies In the data 

nated°so thlt^IGTEr waMbim!!'?*' These J nc0 " ! 'stenc1es were not ellml- 
ZV r DI ? T r” s capability for generating a steady-state balance could 
be tested. For both operating points the converged data have nearly zero 
errors with the Iteration (state) variables del nj adjusted accSrdlnJlJ 

Data for the wet design point are shown In flaure 14fd1 Here »*a4n +ho 
errors ore close to zero. OIGTEN Iterates, to on oper^lno pilnfthat ll clSs. 
to the dry design operating point. Figure 14(e) shows data for the wet off 
design operating points. The Input data do nit give as gSd aS engine ballnce 

SisJgi $o*t£ 90 P ° 1nt f ° r th ® $ame reason as d1scu " sed for^he dryoff - 6 

. . F1 f ure shows the correction coefficients calculated by DIGTEM for the 
dry design point, since they all are close to 1.0, the modeling eauatliSs 
fairly accurately describe the design point. In fact one or more lira? 1st 
ferances from 1.0 would Indicate either bad ln t of one "Le 
Inadequate component models. r more 
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APPENDIX F 

TURBOSHAFT ENGINE MODEL 


DIGTEM Is generalized In the aerothermodynamlc treatment of components. 

It also trims calculations to match a design point* These features can make 
It a useful tool for developing simulations of specific engines having the 
same two-spool, two-stream configuration. Also variations of the turbofan 
engine configuration such as a turbojet or turboshaft can be simulated with 
minor modifications to the Fortran coding. With more extensive modifications 
to the coding, arbitrary configurations can be modeled. 

To demonstrate this capability, a turboshaft engine model was Implemented 
by using DIGTEM. A computational flow diagram of the engine Is shown In 
figure 16. Comparing figure 16 with the two-spool, two-stream engine computa- 
tional flow diagram of DIGTEM In figure 2 Indicates the need to make the 
following changes to the basic DIGTEM model: 


(1) The Inlet model must be eliminated. 

(2) The fan must be eliminated. 

(3) The duct must be eliminated. 

(4) The low-pressure-turbine cooling bleed must be eliminated. 

(5) The low-pressure turbine (l.e., power turbine In. the turboshaft) must 
be disconnected from the fan and connected to a load. 

(6) The nozzle must be eliminated. 

(7) The back pressure on the power turbine must be fixed (at atmospheric 
pressure) with turbine flow (and energy) dumped to atmosphere* 

The turboshaft engine model was Implemented In DIGTEM by 

(1) Using the normalized component maps already In DIGTEM 

(2) Specifying a new design point with Input data -satisfying the following 
conditions: 


P 


2.2 


st 




T 


2.2 



T 


2 


w 


13 


0 





W 6 - w 4.1 


T 


7 


T 


6 


n AB a 


0 


A 0 . A e = 0 
FVGP * CVGP - 0 


(FI) 

(F2) 

(F3) 

(F4) 

(F5) 

(F6) 

(F7) 

(F8) 
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(F9) 


W BLLT " 


0 


(3) Modifying the coding In OIGTEM as follows: For the turboshaft the 

state variables and derivatives are 


VS(1) - XNL 
VS< 2) » XNH 
VS(3) - M3 
VS(4) » T3 
VS( 5) • W4 
VS(6) . 14 
VS( 7) . W41 
VS(8) . T41 


V00T(1) - DXNL 
VD0T(2) . OXNH 
VD0T(3) = DM3 
VD0T(4) « 0T3 
VDOT(5) = DW4 
VOOT( 6) « DT4 
VDOT(7) = 0W41 
VD0T(8) . DT41 


These are the first eight state variables In the state variable list for the 
turbofan engine, and thus no recoding Is needed to set up the state vector and 
the state derivative vector. By setting N • 8 In the main program OIGTEM, 

the order of the system Is specified and the 8x8 Jacobian error matrix-wlll be 
generated. 


.... Recoding of OIGTEM routines was required to account for the aforementioned 
differences In the configurations. A fixed value of load torque was 

set In DSGNPT and was sized to zero the rotor speed derivative at the design 
point. Also correction coefficients CC( 1 4) and CC( 1 6) were redefined In DSGNPT 
to reflect the changed coding In the engine routines. The numerator of CC( 16) 
(eq. (8141)) was set to the load torque. Some coding had to be added In the 
power turbine discharge. That Is, temperature T& was calculated Implicitly 
from the calculated turbine discharge enthalpy h 6 . Convergence was obtained 
by guessing T 6 , using T 6 to compute the h 6 through PROCOM, and then com- 
paring h 6 with calculated h 6 . CC(14) was used to Insure a match at the 

design point. 


Finally recoding was done In subroutine TMRSP, where the Inputs to the 
model were specified as functions of time (open-loop control). For the 

turboshaft engine the Inputs are fuel flow w fA to the main combustor and load 

torque Q loa d change on the power turbine. TMfiSP was set up to give a steo 
change In both fuel flow and load torque. 


Figure 17 shows the transient response of the turboshaft engine to simul- 
taneous steps In fuel flow and load torque. Shown are normalized values of 

fuel flow Wp load torque Qi oa( j> low rotor speed N^, high rotor speed 

N h , combustor pressure P 3 , and turbine Inlet temperature T 4 . Note that Nu 
P3, and T4 all Increase with the addition of fuel. Normally Ni would 
Increase also, but the Increase In load caused Nl to drop off For this 
2-sec transient the Integration time step was 0.01 sec. The printout Interval 
was 0.1 sec. The CPU time was 1.06 sec on the IBM 370/3033 computer. 


Thus It Is possible to use OIGTEM to model engines other than a two-spool, 
two-stream turbofan engine. The resultant engine model will have a realistic 
aerothermodynamlc treatment of Its components and will be scaled to a user- 
specified design point. 
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TABLE !, - TRANSIENT SPECIFICATIONS IN OJGTEH 


ORIGINAL PAGE IK 
OF POOR QUALITY 





Parameter 

Setting 

function 


Default 

Allowable 

rango 


NQPER 

3 

1 - NTOTAL 

Desired Initial 
operating point 

H 

0.01 

>0.0 

Integration time stop, 
soc 

TMAX 

20. 

0.0 for steady 
state, >0.0 fop 
transient 

Ooslrod translont 
length, sec 

TOUT 

Or4 

>H 

Printout Interval, 
sec 

IBDINT - 

1 

0 for explicit, 

1 for Implicit 

Integration method 
selector 

IHPCNV 

0 

0 to use logic 
to generate a new 
matrix, 1 to 
generate a new 
matrix every time 
point 

Matrix update selector 

N 

16 

>0 

System order 


•i 


TABLE II. - BACKWARD- DIFFERENCE INTEGRATION SETTINGS 


Parameter 

Setting 

Function 


Default 

A1 lowable 
range 


V0ELTA- » 

0.001 

>0.0 

Initial perturbation 
of guesses, percent/100 

FRAC 

0.25 

>0 

External control of 
Iteration step size 

ton 

0.001 

>0 

Bottom limit on error 
tolerance for matrix 
linearity, percent/100 

T0L2 

0.01 

>TOL1 

Top limit on error 
tolerances for matrix 
linearity, percent/100 

TOISS 

0.0005 

>0.0 

Solution tolerance, 
percent/100 

MPAS 

50 

>0 

Maximum allowable 
Iteration passes 

T0LPCG 

0.5 

>0.0 

Switch for calculating 
a new matrix 

NOBUG 

0 

0 for 

no debug, 1 
for debug) 

Debug selector 
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Station: 


ORIGINAL FAGu \Q 
OF POOR QUALI7Y 


High-pressure 



Figure L - Schematic of augmented turbofan engine. 


































ORIGINAL 
OF POOR 


PAGE m 
QUALITY 


T^f CP J 

foro.oi 


±$UM- 


0,0851 
0.0000 
fO.OOOO 


0.3OO0 o.onoo 
0.0000 0,0000 
10,0000 10.0000 
84,6000 26,0 
0,01/3 0.0)00 

0,00|) 0,0000 
If). 0000 10,0000 
84.000(1 8!l,ll 
0.01/;* o.oioo 
0.001) 0,0000 
10.0000 10,0000 
04.01100 84,0 
0,08/0 
(1,0080 
1 !) . 0000 

04.0000 H.0 

0.0.WO 0.0, HO 

0.0030 9,0000 

If), 9000 14.0000 

04.0000 89.0 

0,044,’* 0.0390 

0,0044 0.0000 

H.OOOO 14.0000 
,'4.0000 (*10.0 
0,041)4 0.9411 

0.0044 0,0000 

10.0000 14.0000 

(.‘4 .0000 20.0 

0.04/0 0,0423 

0.004/ 0.0000 

13.0000 14.0000 

24.0000 20.0 

0.0400 0.0432 

0.0040 0.0000 

13.001)0 14.2/411 
23.0005 25.0 
0.0490 0.0432 

41100031 >0.0001 
9.5144 11.0440 
23.4514 25.0 
0.0404 0.0409 

**0.01/3 "0.0242 
2.3/93 4.3/93 

20.3/93 22.3/93 
- 0.03*»2 0.0314 

"0.0373 -0.044/ 
0.0000 1./004 

10.5043 11.4/34 
0.0049 0.0000 

-0.03/2 -0.0429 
0.0000 1.0000 
9.0000 10.0000 
0.0000 -0.0024 
-0.032/ -0.03// 
0.0000 1.0000 
9.0000 10.0000 
0.0000 -0.0023 

-•fl rl r*’**** 


WIMI <«•«.<» 




0.7000 

0,9400 

0.7600 

1.0000 

0,0000 

1,0900 

0,0200 

1,1999 

0.0376 

0.0900 


17,0000 

18.0000 

19.0000 

20,0000 

21.0000 

22,0(100 

23,0000 

0,01*0 

0,0121 

0,0104 

0,0006 

0,0069 

0,0052 

0,0034 

17.0000 

18,0000 

19,0000 

20,0000 

81,0000 

28,0000 

23,0000 

0.0137 

0,0170 

0,0)03 

0,0006 

0,0069 

0,0052 

0.0034 

17 0000 

10,0000 

19,0000 

20,0000 

21,0000 

22,0000 

83,0000 

0,0223 

0,0199 

0,016? 

0,0139 

0,0)12 

0,0004 

9,0056 

17,0000 

10,0000 

19.0000 

20,0000 

81,0000 

82,0000 

23,0000 

0,0302 

0,0269 

0,0227 

0,0109 

0,0191 

0.0113 

0.0076 

17,0000 

18.0000 

19,0000 

20.0000 

21.0000 

22.0000 

83.0000 

0.0334 

0,0309 

0,0269 

0,0221 

0,01/7 
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Figure 4 - Concluded. 
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0.1 

0.0 

Z1 

VALUES - CURVE 1 

0.15 

0.15 

0.10 

0.05 

0.00 

Z2 

VALUES • CURVE 1 

.225 

0.225 

0.150 

0; 075 

0.000 

Z3 

VALUES - CURVE 1 

0.0 

0.4 

0.5 

0.6 

0.7 

X 

VALUES • CURVE 2 

0.6 

0.6 

0.4 

0.2 

0 0 < 

Z1 

VALUES - CURVE 2 

0.30 

0.30 

0.20 

0.10 

0.00 

Z2 

VALUES - CURVE 2 

0.450 

0,450 

0.300 

0.150 

0.000 

Z3 

VALUES - CURVE 2 

0.0 

0.6 

0.7 

0.8 

0.9 

X 

VALUES - CURVE 3 

0.9 

0.9 

0.6 

0.3 

0.0 

Z1 

VALUES ** CURVE 3 

0.45 

0.45 

0.30 

0.15 

0.00 

Z2 

VALUES - CURVE 3 

0.675 

0.675 

0.450 

0.225 

0.000 

Z3 

VALUES - CURVE 3 


NCV - NUMBER OF CURVES IN MAP 

NPT - NUMBER OF POINTS PER CURVE 

NFCT - NUMBER OF COMMON FUNCTIONS OF X, Y 

.JCOM - SWITCH FOR COMMON CURVE BREAKPOINTS 

Figure 5. ** Example of map Input data. 


pm * -m 




0 5 10 IS 


Time, sec 

Figure 6. - Open-loop controls for DIGTEM test case* 
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iih 




;^ 3 » 




« s ys R 2ssi N ^s , ? Rsp< mf<v * wr? * cvop ^ , - 

c * A8#TIME* JSS) 

c 

C 


IF JSS .BQ. 0) RETURN 
IF (TIME .GT, 10.0) GO TO 100 
WFAe , 37+TIME*( 1 . 7- . 37 )/2 , 0 
IF (MFA ,GT. 1-. 7) WF4«1.7 
MF7eQ . 0 

FVGP°-24 . 99+(TIME-l , 0 )*(24 , 99-1 . 7 )/l . 00 
IF 6FVGP .IT. -29.99) FVGP«-24.99 
IF (FVGP .GT. -1.7) FVGPa-1.7 
FVGP»BSFVGP-FVGP 
CVGP«-20 , 0+TIME#(20 . i4 . )/2 . 0 
IF (CVGP .GT. 4,0) CVGP«4,0 
CVGP«BSCVGP-CVGP 
A8 e 430 . 0 
AE S 49? . 0 
RETURN 
100 CONTINUE 
WF4 e l . 7 

MF7«0.O*(TXME-1O.O)#C5.O-O.O)/3.O 

IF (MF7 .GT. 5.0) WF7*5,0 
FVGPa-1 . 7** (TIME-1 0 . 0 )*(2 . 5-1 . 7 )/3 . 0 
IF (FVGP .IT* -2,5) FVGpa-2 .5 
CVGPa4, 0 

CVGPaBSCVGP-CVGP 

FVGPaBSFVGP-FVGP 

A8»430, +( TIME- 10, 0)*( 660 ,-430 , )/3. 0 

IF (A8 .GT. 660.) A8«660. 

AE C 492. *■ (TIME-1 0 . 0)#(880 .-492. )/3. 0 
JFJ (AE .GT. 880.) AE*880. 

RETURN 

END 


Figure 7. - Subroutine TMRSP for test case. 


i.'u.CUiVtf. IfVu.iL’ 

OF POOR QUALITY 


O0QQ1QO 

0000200 

0000300 

QQQQ6QG 

0000500 

0000600 

0000700 

OOOOAOO 

0000900 

0001000 

ooonoo 

000)200 

0001300 

ooonoo 

0001500 

0001660 

ooonoo 

ooonoo 

0001906 

0002000 

oooaioo 

0003200 

0002300 

0002600 

0002900 

0002600 

0002700 

0002000 

0002900 

0003000 

0003100 

0003200 

0003300 

0003600 

0003500 

0003600 

0003700 

0003000 

0003900 

0006000 

0006100 

0006200 

0006300 

0006600 

0006600 

0006600 

0006700 

0006000 

0006700 

0005000 

0005100 

0005200 

0005300 

0005600 

0005500 

0005600 

0005700 

0005800 

0005900 

0006000 

0006100 

0006200 

0006300 

0006600 

0006500 

0006600 

0006700 

0006000 

0006900 

0007000 

0007100 

0007200 

0007300 

0007600 

0007500 

0007600 

0007700 

0007800 

0007900 

0008000 

0008100 

0008200 

0008300 

0008600 

0008500 

0008600 

000870C 

0008800 

0008600 

0009060 

0009100 

0009200 

0009300 

0009600 

0009500 

0009600 

0009700 

0009800 

0009900 

0010000 


MAIN PROGRAM POR D I u T E M 
DOUBLE PRECISION EMAT.PETERM 
REAL KOIWHT ,KDLMLT 

DIMENSION LW( 16) ,MW()6) ,VMAT( 16) ,YYYt 16) ,PELEt 16) ,RR( 16) 

COMMON /HELP/ MATRIX 
COMMON/OUTPT/MATTOT 
CQMMON/'SYTURD-'H.IIIPCNV 

^ COMMON XNMAPS/^H JRR^PjMSER) ,pj(8JA> ,FA(9lfl) ,PS(82A), pjtggA), „ 

« c cc< §0 b’To ct AHC , Dr.T AO : SHt A A o l 3 * vs * ' V41 ' W6 ' V7 * XXH ,xx L * oa^vop. B9CVOP , - 

llSISSSIISi** 

! WS^Sm^samarn ; 

? }#DW06#DW13a OW3aDW6 »DW61aDW6#DW7rDXNH#DXNL*BTA#iB. PTA ft PTAHTM - 

S ma IfKBLWHT » KBIWLT # PE* PQ » P0A»P0QT7 »P13*P2»P2A»P22.P9I>0!>m 6* 

WRI TE( 6 dool) ITRANS ’ H0PER ’ 1B01NT . H, TMAX, TOUT, NOBUO, TIME 
1008 FORMATt 1 HI > 

WRITE(6 , 1003) 

„ WKITEC6 » 1006 ) 

J5S? [JORMAT ( 6 OX » 1 KXNKM NtfMIftt KKtHt D I 0 T E M KKKKNKKMKttltKMMtt * 1 

1006 FORMATt 3 OX, -TURBOEAN ENGINE MODEL’) ’ 

WRiTE<6a502) 

WRITE(6 * 1007) 

1007 EORMATUQX, ’INPUT DATA*) 

502 EORMAT(Z) 

N0PER=3 

IOOlNT-1 

Hs*01 

TMAXS20.0 

TOUT=.l 

H = 16 

XHrCNV=0 

C#*»KKREAD MAP DATA 

RCAD( 5» 506 ) KBLWHT »RBtWLT 
CALL IHDATA(N1 » El ) 

READ (5» 505) BSFVGP, BSCVOP 

RRHAREAD^OPERATINO^POINT DATA UP=l DRY DESIGN, IP=NDRYU WET DESIGN) 

NTOTAL=NDRY+NAUG 
DO 100 1P-1.HOPER 
READ(5,507 ) POINT 

CMWMKKKItICKNKHKMMMNKMMMKKKtfMMMKNMIMItKKMNKMK 

c - 0H VAHE 0E0MEm 

CVOP*BSCVGP-CVGP 

C 

IP tIP .ED. 1) GOTO A» 

IWETsNDRYU 

JP <IP -ED. IWET) GOTO AV 
GOTO 50 
C 

* P2*P13»P22aP3 oP6*P61aP5aP6*P7*TAM»T2»T13a T22* T3 
2 ETAAB 1 EN 4 XlV'xNH , ljF4 3 |jlFJ 2 ?* M J^»?^S2 4l ' W0B ' W07 ’ DH4 ' DH ‘ , ^* r A B ' " 

p ! Simlt rKASw*??": A IfttcS x,w,c,,M • ew,, • , * wo, ** - 


SO CONTINUE 

....STEADY-STATE 

ITRANS=0 

CALL ENGINE 


Figure 8, - Main routine DIGTEM tor the output test case. 



LlRtGlNWb 


0010100 

0010800 

0010300 

001040Q 

0010500 

0010600 

0010700 

0010600 

0010900 

0011000 

0011100 

0011200 

0011300 

0011400 

0011500 

0011600 

0011700 

0011600 

0011900 

0012000 

0012100 

0012200 

0012300 

0012400 

0012500 

0012600 

0012700 

0012800 

0012900 

0013000 

0013100 

0013200 

0013300 

0013400 

0013500 

0013600 

0013700 

0013800 

0013900 

0014000 

0014100 

0014200 

0014300 

0014400 

0014500 

0014600 

0014700 

0014600 

0014900 

0015000 

0015100 

0015200 

0015300 

0015400 

0015500 

0015600 

0015700 

0015800 


IPUP .NC. NOPER) GOTO 100 
MRlTG(6tBl4) POINT 

100 CONMNUC 

KMNHMHMMMKMMHMHMMMMMMMNMKflKMNMN 


NMMNMN 


INITIALIZE STATE AND DERIVATIVE VARIABLES 


VS<1) 

VS<2> 

VS(3) 

VS<4) 

VS(5> 

VS(6) 


XHL 

XNH 

W3 

T3 

W4 

T4 


VS(7) a W41 
V$<8) « T41 
VS<9) « W6 
V S C IP) # T6 
VS<11) 8 W7 
VS(12) 8 17 
VS( 13) a WA13 
V5< 14) « 1406 
VSU5) a W13 
VS( 16 ) 8 T13 


VDOT(l) 

VD0TC2) 

VDOT<3) 

VDOT (4) 

VD0T<5) 

VDOT(6) 

VDOT(7) 

VD0T(6) a 

VDO 1(9) a 

VDOT ( 10 > 

VDOTUl) 

VD0TU2) 

VDQTU3) 

VDQTU4) 

VD0TC15) 

VDCTU6) 


DXNL 
DXNH 
DW3 
DT3 
0W4 
DT4 
DW41 
DT41 
DW6 
* OT6 
D147 
OT7 
DWA13 
DW06 
DW1 3 
OT13 


IF ( IBDINT *EQ, 1) OOTO 2000 
CALI FDINTG 
GO TO 3000 
C 

2000 CONTINUE 

CALL BDINTO 
3000 CONTINUE 

504 FQRMAT(5F12, 5) 

505 FORMATC2F10 ♦ 0) 

507 FORMAT(9X,I3> 

506 FORMAT ( (1X»2(12»2X) ) ) 

510 P0RMAT(6F12* 5) 

511 PORMATt3F12, 5) 

514 FORMAT (/* * OPERATING PO IN-?— NUMBER- 
STOP 
END 


, 14 ) . 


Figure a - Concluded, 


INPUT DATA 

OPERATING POINT NUMBER 3 

rimn » o^oooo seconds 


maimi nar&: ^ 

of pooh qmuvf 



STA 2 

STA 13 

STA 8.8 

fiTA 3 

STA 4 

STA 4.1 

STA 6 

STA 7 

PRESSURE 
TEMPERATURE 
DERIVATIVE 
MASS ri!ow 

14.6960 

19.1000 

^rsoo” 

1,98573^^ 

54.0000 

80.6036 

476.80ft 

49.7557 

99.3999 

966.000 

53.3888 

41.5667 

J4.S999 

1560.00 

-8slil95r nn "” 

J7,0999^ 

-.416868"" 

49.3406 

17.5000 

765.000 
io"n36""""" 

104.000 

17.3000 

765,000 

-93.3186 

104.451 

^DERIVATIVE* 


0*844141 B- OS 





0,866595i;-08 


stored s;;r 


4.54641 


0.466166 

0.871859 

1.51610 

1.61108 

1.46046 

DERIVATIVE 


- . If 8446 


771948E-01 

0.7489996-01 

0.468494E-01 

- . 3112790-01 

-.490962 

ENERGY OCR , 


6.93660 


84.8583 

-6 .48174 

-.631104 

19.0407 

-138.146 

DELTA H 





101.309 

27.1966 




low SPIED SPOOL 8 6175.00 RPM 
DERIVATIVE a ‘85,4913 RPM/SEC 

MAIN COMBUSTOR FUEL FLOW *0. 370000 

BLEED MASS PIOUS- 

LOW PRESSURE *0.688276 
HIGH PRESSURE 8 7.50446 
OVERBOARD *0.113444 


HI OH SPEED SPOOL a 943V. 00 RPM 
DERIVATIVE a 16.1689 RPM/SEC 

AFTERBURNER FUEL FLOW *0.000000 

VARIABLE GEOMETRY — 

PVOP * -84 .9900 
CVOP * -80.0000 
THROAT AREA a 430.000 


FSHIFT *0 . 8I6844E-04 
CSHIFT 8-.339364E-08 


CONVBROED STEADY STATE POINT 


TIME * 0.0000 SECONDS 


— 

— 
STA 8 

STA 13 

STA 2,8 

STA 3 

STA 4 

STA 4,1 

STA 6 

STA 7 

PRESSURE 

14.6960 

19,0867 

80.6006 

99.6367 

94,8361 

87.1478 

17.4849 

17.8840 

TEMPERATURE 

518.670 

571.375 

578.098 

966.401 

1578.93 

1115.64 

784,934 

784.934 

DERIVATIVE 


-.461090E-08 


-.133648 

0.844831 

-.87Q769E-0I 

- . 525365E-0I 

0.221376E-01 

MASS FLOW 

103.927 

54.0319 

49.8948 

41.6306 

48.0007 

49.5815 

104.183 

104.183 

DERIVATIVE 


- . 166016E-01 





-♦ 131836E-0I 


STORED MASS 


4.54076 


0.467103 

0.878181 

1.52035 

1.80961 

1.47928 

DERIVATIVE 


915527E-04 


0 . 305176E-04 

- . 562072E-04 

152588E-04 

- . 152588E-04 

0 . 106812E-03 

ENERGY DER. 


- . 209370E-01 


624846E-01 

0 . 666237E-01 

-.138387 

-. 950707E-01 

0 . 327466E-01 

DELTA H 





101,306 

87.3107 

1 

L 


LOW SPEED SrOOt a 6181.88 RPM 
DERIVATIVE a 0.151316E-01 RPM/SEC 

MAIN COMBUSTOR FUEL FLOW * 9.370000 

BLEED MASS FLOWS— 

LOW PRESSURE *0.689648 
HIGH PRESSURE * 7.58079 
OVERBOARD -0.113691 


HIOH SPEED SPOOL = 9444.64 RPM 
DERIVATIVE = 0 .864081E-01RPM/SEC 

AFTERBURNER FUEL FLOW *0.000000 

VARIABLE GEOMETRY — 

FVGP * -84.9900 
CVGP a -80.0000 
THROAT AREA = 430.000 


FSHIFT «0*8l8S38E-04 
CSHIFT *- .415303E-02 


TIME 

PO 

TAM 

XHt 

P2 

T2 

XNH 

P13 

T13 

WP4 

P22 

T22 

WF7 

P3 

T3 

A8 

P4 

T4 

AE 

P41 

T41 

FVGP 

P6 

T6 

CVGP 

P7 

T7 

MATTOT 

0.000 

14.696 

516.67 

6181.8 

14,696 

516.67 

9444.8 

19.069 

571,37 

0.37000 

80,601 

576.10 

0.00000 

99.637 

966.40 

430.00 

94.836 

1578,9 

498.00 

87,148 

1115.6 

-84,990 

17,465 

764.93 

-80,000 

17.864 

764.93 

1 

0.190 

14,696 

516,6? 

6192.8 

14.696 

518.67 

9463,8 

19,196 

578.50 

0.43650 

80.687 

576,96 

0,00000 

108.31 

973.81 

430.00 

97.634 

1664.0 

498.00 

87.936 

1171.0 

-84.990 

17,630 

805,78 

-16.800 

17.480 

601,58 

8 

0.800 

14,696 

518,67 

6835,3 

14.696 

516,67 

9513.4 

19.359 

574.89 

0.50300 

80.835 

560.88 

0.00000 

106.84 

963.78 

430.00 

101.57 

1786.0 

498.00 

86 , 988 
1845.0 
-84.990 

17.798 

639.66 

-17,600 

17.568 

655.75 

8 

0.300 

14.696 

518.67 

6309.9 

14.696 

516.67 

9597,1 

19.580 

576.45 

0,56950 

81.053 

561.93 

0.00000 

111.44 

997,58 

430.00 

106.67 

1875,5 

498.00 

30.059 

1308,1 

-84.990 

17,969 

869.95 

-16.400 

17,748 

666,61 

3 

0.400 

14,696 

518.67 

6416.8 

14.696 

518.67 

9704.8 

19.893 

579.81 

0.63600 

81,365 

584.41 

0.00000 

117.80 

1014.4 

430.00 

118.66 

1958,8 

498.00 

31.467 

1362.7 

-84.990 

16,868 

696.98 

-15.800 

16.006 

694,26 

3 

0.500 

14.696 

518.67 

6551.6 

14.696 

518.67 

9686.6 

80.896 

582.74 

0.70850 

81.818 

587,76 

0.00000 

185.16 

1038,5 

430.00 

119,97 

8017.4 

498.00 

33.104 

1408.3 

-84.990 

18.610 

980.44 

-14,000 

16.387 

918.49 

3 


Figure 9* -DIGTEM output for test case. 
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ORIGINAL, PAGE IS 
OF. poor QUALITY 


3,300 

14.696 
516,67 
9254 . 3 

14,696 

516.67 

11923 , 

34,138 

698.47 

1,7000 

35,685 

698,42 

0.00000 

269.48 

1327,9 

430,00 

254,57 

25 , 90.7 

492,00 

69,894 

1767,8 

- 1,7000 

31,253 
1167 ,8 
4 ,0000 

30.662 

1157,7 

4 

3,400 

14,696 
81 A . 67 
9288,1 

14,696 

816,67 

11928 . 

34,140 

698,50 

1,7000 

38.687 

698.45 

0.00000 

265.54 

1327.7 

4 . 30.00 

254,62 

2530,6 

492.00 

69,608 

1787,4 

" 1.7000 

fura 

4,0000 

30.668 
1187. 6 
4 

3,800 

14,696 

816,67 

9258,7 

14,696 
516,67 
11927 , 

34,143 

695,53 

1,7000 

38,689 

698,47 

0,00000 

265,56 

1327,8 

430,00 

492 , 4 o 

69,618 

1767,4 

- 1,7000 

l ! 87?6 

4,0000 

30,674 

1187,9 

4 

3,600 

14,696 

9286 ?! 

14,696 
816,67 
11926 , 

34,143 

695,58 

1,7000 

35,689 

696,49 

0,00000 

265,61 

1927,9 

4 , 90,00 

2930,5 

492,00 

69,622 

1767,4 

- 1,7000 

31,288 

1197,9 

4,0000 

30,673 

3,700 

m 

9286,6 

14,696 

816,67 

11929 , 

34,144 

695,87 

1,7000 

39,690 

698,90 

0,00000 

269 , 6.9 

1326.0 

430.00 

254.71 

2930,5 

492,00 

69 , 6*7 

1767,3 

- 1.7000 

31,689 

1196.0 

4.0000 

1M!J 

3 , A 00 

14,696 

816.67 

9266,9 

14,696 

816,67 

11929 . 

34,144 

698,98 

1,7000 

39,690 

696,91 

0,00000 

264,65 

1328,1 

430,00 

nh’s 

492,00 

69.631 

1767,3 

- 1.7000 

31,299 

1196.0 

4.0000 

30,672 

1198,0 

3.900 

14,696 

816,67 

9287,2 

14,696 

516,67 

11930 , 

34.146 

695,60 

1,7000 

35.692 

698,92 

0.00000 

263,66 

1328.1 

430,00 

234,73 

2330.4 

492.00 

69.634 

1767,3 

- 1,7000 

31,294 

1196.0 

4.0000 

30,676 

1196,0 

4 

4,000 

14,696 

516.67 

9267.3 

14,696 

518.67 

11930 , 

34,146 

698,60 

1.7000 

39.692 

696,93 

0,00000 

263,66 

1328.1 

430,00 

234.74 

2930.3 

492.00 

69,636 

1787,3 

- 1,7000 

31,296 

1196.0 

4.0000 

30.677 

1196.0 

4 

4,100 

14,696 

516.67 

9287,4 

14,696 

918.67 

11930 . 

34,146 

699,60 

1,7000 

35,692 
698.93 
0 ,00000 

263,67 

1326.1 

430,00 

234,74 

2330.4 

492.00 

69.637 

1787.3 

- 1.7000 

31,299 

1198.0 

4.0000 

30.676 

1198.0 

4 

4,200 

14,696 

516.67 

9257.6 

14.696 

516.67 

11930 , 

34.146 

695.60 

1,7000 

39,691 

698.53 

0.00000 

263,67 

1326.1 

430,00 

294.74 

2930.6 

492.00 

69.638 

1767,3 

- 1,7000 

31.263 

1138.0 

4.0000 

30,672 

1198.0 

4 

4.300 

14.696 

816.67 

9257.6 

14.696 

518,67 

11930 . 

34,147 

695.61 

1.7000 

39.692 

698.54 

0.00000 

263,68 

1328.2 

430.00 

254,73 

2330.4 

492.00 

69.639 

1787,3 

- 1.7000 

31.299 

1198.0 

4.0000 

30.674 

1198.0 

4 

30.676 

1198.0 

4 

4,400 

14.696 

516.67 

9257.7 

14.696 

918.67 

11930 . 

34.147 

695.62 

1.7000 

39.693 

696.54 

0.00000 

265.68 

1328.2 

430.00 

234.79 

2930.4 

492.00 

69.639 

1767.3 

- 1.7000 

31.298 

1198.0 

4.0000 

4.500 

14.696 

516.67 

9257.7 

14,696 

516.67 

11931 . 

34.147 

695.62 

1.7000 

38.693 

698.54 

O.OOOQO 

263.66 

1328.2 

430.00 

234.79 

2330.4 

492.00 

69.640 

1787.3 

- 1.7000 

31.299 

1198.0 

4.0000 

30.677 

1198.0 

4 

4.600 

14,696 

518.67 

9257.7 

14.696 

518.67 

11931 . 

34.146 

695.61 

1.7000 

39.692 

698.54 

0.00000 

263.68 

1326,2 

430.00 

254.76 

2330.4 

492.00 

69.640 

1767.3 

- 1.7000 

31,263 

1158.0 

4.0000 

30.674 

1198.0 

4 

4.700 

* 4 .696 
518.67 
9257,7 

14.696 

518.67 

11931 . 

34.147 

695.61 

1.7000 

35.692 

698.54 

0.00000 

265.68 

1328.2 

430.00 

254.76 

2330.4 

492.00 

69.640 

1787.3 

- 1.7000 

31.262 

1138.0 

4.0000 

30.674 

1138.0 

4 

4 .800 

14.696 
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6 

13.600 

14.696 

310,67 

9223.5 

14,696 

518.67 

11909. 

34.411 

696.12 

1.7000 

35.920 

698.37 

5,0000 

266.63 

1326.2 

660.00 

255.65 

2523.2 

660.00 

69,696 

1762,3 

-2.5000 

31,631 
1160. 6 
4.0000 

31.056 

2679.1 

6 

13.700 

14.696 
S\8 . 6? 
9223,3 

14.696 

318.67 

11909, 

34.411 

696.1? 

1.7000 

35.920 

698.3? 

5.0000 

266.63 

1326.2 

660.00 

255.65 

2523.2 

600.00 

69.696 

1782.3 

“2.5000 

31.631 

1160.6 

4.0000 

31,059 

8679,1 

6 

13.600 

14.696 

510.67 

9223,3 

14.696 

518.67 

11909, 

34.411 

696.11 

1.7000 

35,920 

698.37 

5.0000 

266.63 

1326.2 

660.00 

255.66 

2523,2 

600,00 

69,696 

1782.3 

“2.5000 

31,632 

1160.6 

4,0000 

31.056 

2679.1 

o 

13.900 

14.696 

510.67 

9223.5 

14.646 

518.67 

11909, 

34.410 

696.11 

1.7000 

35.920 

698.37 

5,0000 

266.63 

1326.2 

660.00 

255,66 

2523,2 

660.00 

69,696 
1762.3 
“2 5000 

31,633 

1160.6 

4.0000 

31,056 

2679,3 

6 

16.000 

14.696 

510.67 

9223.5 

14.696 
518.67 
11909 . 

34.411 

696.11 

1.7000 

35,920 

698.3? 

5,0000 

266.64 

1326.2 

660.00 

255.65 

2523,2 

680.00 

69.696 

1782.3 

“2.5000 

31.632 

1160,6 

4.0000 

31.057 

2679,2 

6 

16.100 

14.696 

310.67 

9223.5 

14,696 

510.67 

11909. 

34,411 

696.12 

1.700C 

35,920 

698,37 

5,0000 

266.63 

1326,2 

660.00 

255.66 

2523,2 

660.00 

69.896 

1782,3 

“2,5000 

31.631 

1160.6 

4,0000 

31,058 

2679.1 

6 

16,200 

14.696 

510.67 

9223.3 

14,696 

310,67 

11909, 

34.411 

696.12 

1,7000 

35,921 

698,37 

5.0000 

266,63 

1326,2 

660.00 

255,65 

2523.2 

660,00 

69,896 

1762,3 

“2,5000 

31,630 

1160.6 

4,0000 

31.060 

2679.0 

6 

16.300 

14.696 

318,67 

922 3.5 

14.696 

318.67 

11902. 

34,411 
696.12 
J . 7 n ft o 

35.920 
698,37 
5 , noon 

266.63 

1326.2 

6M,ftn 

255,65 
8523,2 
88 '1,0 0 

69,696 

1762.3 

-2.4PP0 

31,632 

1160.6 

•’.0000 

31,059 

2679,1 

r, 

16.400 

14.696 

318,67 

9223.6 

14,696 

318.67 

11909. 

34,410 
696, a 
1,7000 

35.920 

696.3? 

5.0000 

266,64 

1326,2 

660.00 

255,65 

2523.2 

680.00 

69,696 

1762,3 

-2,5000 

31,635 

1160,6 

4.0000 

31.055 

2679,3 

fi 

16,300 

14.696 

510,67 

9223,5 

14,696 

510,67 

11909. 

34.410 

696.11 

1.7000 

35.920 

698.37 

5,0000 

266.64 

1326,2 

660,00 

255.66 

2523.2 

660.00 

69.896 

1768.3 

“2,5000 

31,633 

1160.6 

4.0000 

31,056 

2679,2 

6 

16.600 

14,696 

518.67 

9223.3 

14,696 

510,67 

11909. 

34.411 

696,12 

1.7000 

35.920 

698.37 

5,0000 

266.64 

1326.2 

66A.00 

255.65 

2523.2 

680,00 

69,696 

1782.3 

“2,5000 

31,631 

1160.6 

4.0000 

31,057 

2679.2 

6 

16.700 

14.696 

518.67 

9223.3 

14.696 

510,67 

11909. 

34.411 

696.12 

1.7000 

35,921 

698.38 

5,0000 

266.64 

1326.2 

660.00 

255.65 

2523.3 

860.00 

69,696 

1762.3 

-2.5000 

31,628 

1160.6 

4.0000 

31,061 

2676,9 

6 





Figure 9. - 
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ORIGINAL PAGE 19 
OF POOR QUALITY 



fH »,67 

W 3.6 

16609 

* 16.67 

11909 . 

n,m 
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11000 , 

U '000 

«*2! 
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moo , 
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W 1 .& 
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11000 , 

17 , m 
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11900 . 

17,500 
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17 ,<100 

19,696 

916,67 

9223.9 
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5 ) 0.67 
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17,500 

19,696 

916,67 

9223,9 

19,696 

910,67 

11909 . 

17,600 

19,696 

916,67 

9223,9 

19,696 

510.67 

11909 . 

17,700 

19,696 

916,67 

9223.5 

19,696 

516,67 

11909 . 

17,000 

19,696 

916.67 

9223.9 

10,696 

516,67 

11909 . 

17.900 

16.696 

916.67 

9223,9 

19,696 

516,67 
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IS . 000 

19.696 
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19.696 

516.67 
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516.67 

9223.9 

19.696 

516.67 
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516.67 
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16.600 
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516,67 

11909 , 

16,900 
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19.696 

916 . 6 ? 

11909 . 
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516.67 
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516,67 
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19.696 
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11909 . 


36.611 
606.12 
1 . 7000 

36.921 

600 . 3 / 

5,0000 

266.63 

1326.2 

660.00 

256.66 

2523.2 

000.00 

69.097 

1 / 02.3 

“ 2.6000 

31,131 

1160,6 

6.0000 

31.060 

2679.0 

6 

36.011 

606.11 
1,7000 

36.920 

690,37 

6.0000 

266,66 

1326.2 

660.00 

256.66 

MM 

60,097 

1702,3 

“ 2.6000 

31,632 

1160.5 

6,0000 

31,069 

2679,1 

6 

36,611 

606,11 

1,7000 

36,920 

690,37 

6,0000 

266,63 

MM 
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000 00 

69,096 
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* 2 , 500(1 

31,632 

11611,6 

0,0000 

31,060 

2679,1 

6 

J 0 , 6 U 

606,11 

1,7000 

36,920 

690,37 

5,0000 

266,66 

1126,2 

660,00 

266,66 

2623,2 

000,00 

69,096 

1 / 02,3 

•* 2.6000 

J ) ,632 
1160,6 
6,0000 

31,060 

mi 

1 7000 

35,920 

690,37 

5,0000 

266,66 

1326,2 

660,00 

266,66 

2623,2 
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69,096 

1702,3 

*' 2,6000 
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6,0000 

M 

6 

36,011 

606,12 

1,7000 

36,920 

690,17 

6,0000 

266 , 6,1 

H 26,2 

660,00 

km 

000,00 

- 2,5000 

31,631 

1160,6 

6,0000 

31,050 

2679,1 

6 

39,611 

696,12 

1,7000 

35,921 

690,37 

5,0000 

266,63 

1326,2 

660,00 

255,66 

2523,2 

600,00 

69,096 

1702,3 

» 2 , 60 tl 0 

31,631 

1 ) 60,6 

6,0000 

31,060 

2679,0 

6 

30,911 

606,11 

1,7000 

36,920 

690 . 3 / 

5,0000 

266,63 

1 , 126,2 

660,00 

266,65 

2623,2 

600,00 

" 2,5000 

31,632 

1160,6 

6,0000 

31,059 

2679,1 

6 

39,611 

696.11 

1.7000 

35,920 

090.37 

6,0000 

266,63 

1326,2 

660,00 

265.66 

2523.2 

600.00 

69,096 

1762.3 

" 2.5000 

31,632 

1160.6 

6,0000 

31,050 

2679.1 

6 

39,911 

696,11 

1.7000 

36,920 

690,37 

6.0000 

266,65 

1326,2 

660.00 

255,66 

2523.2 

600.00 

69,096 

1702,3 

" 2,5000 

31,632 

1160,6 

6.0000 

31,037 

2679.2 

6 

39,911 

696,11 

1,7000 

35.920 

696 . 3 ? 

6.0000 

266,69 

1326.2 

660,00 

255 65 
2523.2 
600.00 

69,097 

1762.3 

" 2,5000 

31.631 

1160,6 

6,0000 

31.050 

2679,1 

6 


36 . 6 U 

696.12 

1.7000 

35.920 
690.37 
5 . mi 6 0 

266.63 

1326,2 

660,00 

253,66 

2325,2 

680,00 

69.096 
1702 3 
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31,631 

1160,6 
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31,058 
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6 
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696.12 

1.7000 

35.920 

696.37 

5.0000 

266,63 

1326.2 

660.00 

258.66 

2523.2 

680,00 

69,896 

1782.3 

“ 2.5000 

31.631 

1160.6 

6.0000 

31,059 

2679.0 

6 

36.611 

696.11 

1.7000 

35.920 

696.37 

5.0000 

266.69 

1326,2 

660.00 

235,63 

2523.2 

680.00 

69.696 

1762.3 

" 2.5000 

31.632 

1160.6 

6.0000 

31.059 

2679.1 

6 

36.611 

696.11 

1.7000 

35.920 

696.37 

5.0000 

266,63 

1326.2 

660.00 

255,66 

2523.2 

660.00 

69.096 

1782.3 

“ 2.3000 

31.632 
1160,6 
6 , 0000 

31.058 

2679.1 

6 

36.611 

696.11 

1.7000 

35.920 

696 . 3 ? 

5.0000 

266.63 

1326.2 

660.00 

253.65 

2523.2 

880.00 

69.696 

1702.3 

“ 2.5000 

31.632 

1160.6 

6.0000 

31.058 

2679.1 

6 

36.611 

696.12 

1.7000 

35.920 

696 . 3 ? 

5.0000 

266.69 

1326.2 

660.00 

253.63 

2523.2 

600.00 

69.896 

1702.3 

“ 2.3000 

31.631 

1160.6 

6.0000 

31.058 

2679.1 

6 
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696.12 

1.7000 

35,920 

696 . 3 ? 

5.0000 

266.63 

1326,2 

660.00 

253.66 

2523,2 

880,00 

69.696 

1762.3 

" 2.5000 

31,631 

1160,6 

6,0000 

31,056 

2679.1 

6 

36.611 

696,12 

1.7000 

33,920 

690.37 

5.0000 

266.03 

1326.2 

660.00 

255,66 

2523,2 

680.00 

69,896 

1782,3 

“ 2.3000 

31.631 

1160,6 

6,0000 

31,059 

2679.1 

6 

36.611 

696.11 

1.7000 

35,920 

696 . 3 ? 

3.0000 

266.69 

1326,2 

660.00 

255.65 

2523,2 

000.00 

69,896 

1702,3 

" 2,5000 

31.638 

1160,6 

6.0000 

31.050 

2679.1 

6 

36.611 

696.11 

1.7000 

33.920 

690 , 3 ? 

3,0000 

266,63 

1326.2 

660.00 

255.66 

2523.2 

880.00 

69,096 

1702,3 

" 2.3000 

31.632 

1160.6 

6,0000 

31.056 

2679,1 

6 
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696.11 

1,7000 

33,920 

690 . 3 ? 

3.0000 

266.63 

1326,2 

660.00 

253.63 
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69.096 

1762,3 
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6.0000 
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6 

36.611 

696,12 

1.7000 

33,920 

690 . 3 ? 

3.0000 

266.69 

1326,2 
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660.00 

235.63 

2523,2 

880.00 

69.896 
1 782,3 
" 2.5000 

31.631 

1160,6 

6.0000 

31,059 

2679.1 

6 

36,611 

696.12 

1,7000 

33.920 

696 , 3 " 

3.0000 

266,63 

1326.2 

660,00 

255.63 

2523.2 

680,00 

69.896 

1782.3 

" 2.5000 

31.632 

1160.6 

6.0000 

31.056 

2679.1 

6 


Hgure 9* 

- Continued. 






ORIGINAL TAki; O 
OF POOR QUALITY 






Combustor pressure. Turbine inlet and augmentor 

psi temperatures. T^and T 7 . °R 






Integration time step, see 

(b) Number of Jacobian matrices. 

Figure 11, - Integration time step study for DIGTEM 
sec test case transient 






Combustor pressure, P3. N/cm 2 low rotor speed. N L . rpm 
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(b) Combustor pressure. 


Figure 12. - Comparison of low rotor speed and combustor 
with different Integration time steps. 


pressure responses for the 


DICTEM test case 
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INPUT DATA 

OPERATING POINT NUMBER l 


0,0000 SECONDS 


_PRmUR0 

jeiiPCRAtufir 
derivative^ 
^JIAOS Flow*" 
^DERIVATIVE* 
^STORED MASS 
DERIVATIVE** 
JNEROY DER?* 

JimT** 


518*670* 


193.000 


OTA 13 

5TA 2,2 

30,5000 

696.000 

-.268700E-01 

36,0000 

697.999 

86,5000 
0 . 390625E-01 

107.000 

6.73723 


-.213623E-03 


-.181032 





^267,000 
^1325*00* 
-.080204 
88. 0990* 


236.000 

jopiToo" 

0.990337* 

89.-7999** 


70.0000 
^1780 , 00**** * 
0;193770C-0r 

107.000 


-.002098 


0.460206 

^3l7970E*0r 

0,255592 

167.000 


J.05708 
;.057760E-00 
0.0761958*01* 
75.5001 


31,8000 

nilTr 

100589"* 

^ 190 . 900 "" 

0^2730 38B- 
2.22702 " 
-,.1220706* 
-.322003** 


JO, 5999 
J160.00* 
-.903200* 
190.902* 


^ 1J7208 

lh 

-.891712 


aryut ** v/dd . nn adm 

DERIVATIVE a 0.722958E-01 RPM/SEC 


MAIN COMBUSTOR FUEL FLOW a 1.70000 


HI0H J£ E6D 3P00L c 11900.0 RPM 
DERIVATIVE a 0 , 309328E-01RPM/SEC 


BLEED MASS FLOWS— 

LOW PRESSURE # 1,03999 
HJ OH PRESSURE a I7,g000 
OVERBOARD »0. 260009 


AFTERBURNER FUEL FLOW =0*000000 


VARIABLE GEOMETRY — 

FVGP a **1.70000 
CVOP a 0.00000 
THROAT AREA a 930,000 


FSHIFT s 0 .398606E-Q6 
CSH1FT *0 .211876E-07 


CONVERGED STEADY STATE POINT 
T *ME a 0,0000 SECONDS 


PRESSURE 


^DERIVATIVE** 
MASS FLOW** 
^ DERIVATIVE* 
STORED MASS* 
DERIVATIVE*" 
^ENEROY*DER?" 
DELTA H 


STA 2 

STA 13 

STA 2.2 

10.6960 

£18.670 

30.0997 

695.999 

-.620850E-03 

35,9997 

698.000 

193.500 

86.5008 
0 . 117188E**01 

106.999 


6.73719 



-.msseE-oo 



-.020976E-02 



266.999 

1325.00 

-.708905* 

88.0983 


255.999 

2520.01* 

1.53016 

89.7990 


69.9996 

1780.00* 


106.999 


0.060202 


2.05706 


0.700206 

167.001 


0.136030 

75.5002* 


v ". y t-tw srugi - TeUU.a!) RPM 

DERIVATIVE a -.361O77E-01 RPM/SEC 


MAIN COMBUSTOR FUEL FLOW a 1.70000 


HIGH SPEED SPOOL a U900.0 RPM 
DERIVATIVE a 0 . 089058E-01RPM/SEC 


STA 6 

STA 7 

31,7996 

30.5996 

1160.00 
0 . 199381E-01 

1160.00 

-.3606016-02 

190.900 

190.900 

-.251065E-01 


2.22700 

-.9155276-00 

1.77206 
-. 152588E-00 

0 . 0000216-01 

-.6O6167E-02 


BLEED MASS FLOWS— 

LOW PRESSURE a 1,63998 
HIGH PRESSURE a 17.1999 
OVERBOARD =0.260008 


AFTERBURNER FUEL FLOW =0.000000 


VARIABLE GEOMETRY — 

FVGP a **1.70000 
CVOP a 6.00000 
THROAT AREA a 030.000 


FSHIFT 8-.66792XE-06 
CSHIFT a 0 . 000000 


ta) Operating point 1 (dry design point). 
Figure 14. - Steady-state operating points. 
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INPUT DATA 

OPERATING POINT NUMBER ? 
TIME » 0.0000 SECONDS 


— — 

1' " 

ST A 2 

13 

STA 8.2 

STA 3 

STA 4 

STA 4,1 

STA 6 

OTA 7 

PRESSURE 

14.6960 

24.2000 

26.1800 

158,700 

191,300 

41.2399 

21.7000 

21,2100 

TEMPERATURE 

318.670 

614.000 

620,941 

1106,00 

1918.00 

1343,90 

891.000 

891*000 

DERIVATIVE 


-12 8615 


8.30729 

18,0497 

30,6131 

10.6789' 

186.017 

M\SS PLOW 

147,379 

75.0000 

72.4762 

60,1436 

60.9108 

78,1240 

148,000 

147.412 

DERIVATIVE 


0,6230002-01- 





-.8080276-01 


STORED MASS 


5,33696 


0.690088 

0.337861 

1.91891 

1.97831 

1.39913 

DERIVATIVE 


970764B-01 


O.289764E-01 

-.2718232-01 

181122E-01 

o.omooe-oii 

0.988333 

ENERGY DER. 


-69.0056 


3.40044 

6.45929 

58,7438 

21.1288 

201.917 

DELTA H 





127.591 

46.7860 




LOW SPEED SPOOL « 7706.00 RPM 
DERIVATIVE # -19.9786 RPM/SEC 

MAIN COMBUSTOR FUEL FLOW «0. 788000 

BLEED MASS PLOWS— 

LOW PRESSURE 60.937253 
HIGH PRESSURE # U.I931 
OVERBOARD = 0. 169234 


HIOH SPEED SPOOL 
DERIVATIVE a 


* 10434.0 
6.62954 


RPM 

RPM/SEC 


AFTERBURNER FUEt FLOW #0. 000000 

VARIABLE GEOMETRY — 

PVOP -24.9900 
CVGP 8 -4.80000 
THROAT AREA - 430,000 


PSHIPT 80.470579B-04 
CSHIPT 8-.141049E-03 


CONVERGED STEADY STATE POINT 
TIME s 0.0000 SECONDS 



STA 2 

STA 13 

STA 2.2 

STA 3 

STA 4 

STA 4,1 

STA 6 

STA 7 

PRESSURE 

14.6960 

24.1867 

26.1354 

198.752 

131.542 

41.2461 

21.7204 

21.2331 

TEMPERATURE 

318.670 

613.302 

620.708 

1106.18 

1917.88 

1343,57 

892.094 

892.054 

DERIVATIVE 


0.418344E-02 


0,175164 

0 . 248939E-01 

-♦ 321751E-01 

-.403699E-01 

0.000000 

MASS FLOW 

146.998 

74.5027 

72.4956 

60.1913 

60.9291 

72.1270 

147.567 

147.567 

DERIVATIVE 


0 . 351563E-01 





-♦388184E-01 


STORED MASS 


3.35836 


0.650)95 

0.357983 

1.91839 

1.97802 

1.39898 

DERIVATIVE 


0 . 915527E-04 


-.274658E-03 

0.222266E-03 

-.457764E-04 

0.305176E-04 

0.000000 

ENERGY DER. 


0.224164E-01 


0.113891 

0.891161E-02 

-.617243E-01 

-.798326E-01 

0.000000 

DELTA H 





127.600 

46,7440 




LOW SPEED SPOOL * 7694.30 RPM 
DERIVATIVE * 0.229609E-0I RPM/SEC 

MAIN COMBUSTOR FUEL FLOW #0. 738000 

BLEED MASS PLOWS— 

LOW PRESSURE «0. 937482 
HIGH PRESSURE * 11.1978 
OVERBOARD 80,169275 


HIGH SPEED SPOOL * 10437,4 RPM 
DERIVATIVE » 0 .796561E-02RPM/SEC 

AFTERBURNER FUEL FUOW *0,000000 

VARIABLE GEOMETRY — 

FVGP a -24.9900 
CVGP a -4.80000 
THROAT AREA # 430.000 


FSHIFT 80.469834E-04 
CSHIPT 8-.296314E-03 


(b) Operating point 2. 
Figure 14 • Continued. 



INPUT PATA 

OPERATING POINT NUMBER 


ORIGINAL PAGE 19 
OF POOR QUALITY 


0.0000 ftf?CONPft 


PRESSURE 


TEMPERATURE I 518,670 


DERIVATIVE 


MASS PLOW | 105.967 


DERIVATIVE 


STORED MASS 


DERIVATIVE 


ENERGY DER. 


DELTA H 


LOW SPEED SPOOL * 6175.00 
DERIVATIVE c -25.6915 


1.52573 


56.0000 


0 . 244141E-03 


4.56661 


17.5000 



RPM 

RPM/SEC 


MAIN COMBUSTOR FUEL FLOW =0.570000 

BLEED MASS PLOWS— 

LOW PRESSURE =0.628276 
HIGH PRESSURE = 7.50448 
OVERBOARD =0.113444 

CONVERGED STEADY STATE POINT 


0.271259 


1.51610 


0.466186 


*. 771 942E-01 1 0 . 742999E-01 | 0 , 462494E-01 


24.8585 -6.48174 


K 631104 


101.309 


HIGH SPEED SPOOL = 9439.00 RPM 
DERIVATIVE # 16.1629 RPM/SEC 

AFTERBURNER FUEL FLOW =0.000000 

VARIABLE GEOMETRY — 

FVOP = -24.9900 
CVGP = -20.0000 
THROAT AREA « 430.000 


1.48046 


•311279E-01 |-. 450562 


*138.146 


FSHIFT =0 . 216844E-04 
CSH1FT 8-.339364E-02 


0.0000 SECONDS 



PRESSURE 14.6960 


TEMPERATURE 518.670 


DERIVATIVE 


MASS FLOW 


DERIVATIVE 


STORED MASS 


DERIVATIVE 


ENERGY DER. 


DELTA H 


- . 166016E-01 


-.915527E-04 


-.209370E-01 


LOW SPEED SPOOL * 6181.22 RPM 
DERIVATIVE = 0.1513I8E-0I RPM/SEC 

MAIN COMBUSTOR FUEL FLOW =0,370000 

BLEED MASS PLOWS— 

LOW PRESSURE =0.629642 
HIGH PRESSURE = 7.52079 
OVERBOARD =0,113691 


0.467103 10.272121 I 1.52035 


0.305176E-04 -.562072E-04 -.152588E-04 


-.624246E-01 | 0 .666237E-01 -.132387 


1.80961 


1.47922 


-. 152588E-04 0.106812E-03 


-.950707E-0I | 0.327466E-01 


101.306 


27.3107 


HIGH SPEED SPOOL = 9444.84 RPM 
DERIVATIVE = 0 . 264081E-01RPM/SEC 

AFTERBURNER FUEL FLOW =0.000000 

VARIABLE GEOMETRY — 

FVGP = -24.9900 
CVGP = -20.0000 
THROAT AREA = 430.000 

(0 Operating point 3. 

Figure 14, - Continual. 


MHIFT B 0.21«S33e-04 
CSHIFT 8-.4I5303E-02 




































































































































INPUT OATA 

OPERATING POINT NUMBER A 


original, page la 

° F POOR qwurp. 


TIME » 


&40&0 9GC0ND9 


^^PRnOfiURE 

JWERATURt! 

OKRtvAnvn" 

^MABS PlOW** 
OCRIVATIVE* 

jroRco^MAar 

derivative"" 

cneroy*dcr? m 

mnT" 


STA 2 


J44969 

91847Q* 


195*900- - 


»TA 13 


54.0000 

t*"n*.or,» wnBBB 

RH.O^O 


l^ijojiie-oi 

ofi T sqoo 

• .miese-ir 

-.ioeiii 


!!'!.' SPEED SPOOL * 9200.00 PPM 

UtRIVAT, 2 « O.JtlO?*e"oi RPM/SeC 

MAIN COMBUS1 OR FUEL FLOW « 1.20000 

IMPED MASS FlOWS— 

LOU PRESSURE a 1.45999 
Ml OH PRESSURE « iyfloOO 
OVERBOARD ^ 0 * 26 0009 


OTA 8 ,ft 


364000 

690.000" 


107.000 


OTA 5 


867.000 

"i!!!***"*' 

00 . 0999 "-" 


0 t 918947 

oJjoJsJe-os 
. moil"" — 


OTA A 


896.000 

898040" 

ojoojxr 

osrjjio" 


0.460846 
;.568Me-05 
i?267I42 11 

167.000 


OTA 44 


704006 
176040 
0407796041 

IO 7 .OO 0 ""*" 


249740 ^ 
410598E-04 
04677655-01* 
79,9001 


CONVERGED STEADY STATE-POINT 
1 IMK 2 04000 SECONDS 


H *OW SPEED SPOOL 0 11900 4 ppm 
DERIVATIVE * 0 .698657E-02RPM/9EC 

AFTERBURNER FUEL HOW * 940000 
VARIABLE GEOMETRY — 

JVOF * -8,90040 
CVOP # 440000 
THROAT ARlUJII 6*0.000 


OTA 6 


314000 
116040**"" 
;4449B9*"*" 
J94440 
*4125006-01 
248708*""" 
;428070 e"o 5 
.388005 


OTA 7 


J04000 

**866840 


* 475999n*»0i 
196,009*"""' 


0406198 
145141 M 
. 459493E41 


FONT FT *-4388222-02 
CSHIPT *04116761-07 


PRESSURE 

ftMPtRATURE 


I DERIVATIVE 


j m/.ss now 
I dertvative" 


MASS 


DERIVATIVE 
tNCROr DER. 

deita*h""* 


STA 2 

STA 13 

STA 24 

STA 3 

14.6960 

518.670 

34.6747 

696497 

*4597106-02 

34.1188 

697491 

267469 

1323.39 

*4537466**01 

68.3327 

193.139 

85.8615 
*♦ 390625E-02 

107478 





6.76653 

0406812E-03 

*.311064E-01 


0.915664 ( 
0.12310SE-0J • 
-.87S3UE-01 C 


^DERivATIvn^roI^^OlP Rp8&C 

MAIN COMBUSTOR FUEl FLOW » 140000 



OTA 44 


704325 

177741 

430177E-01 

107474 


846686 

^762939E-04 

159419 

794431 


BLEED MASS FLOWS— 

LOW PRESSURE * 144340 
HI ON PRESSURE « 17.2407 
OVERBOARD *0 460629 


MIOM Blp.u*r?l!2 0L * IMMMM rpm 
DERIVATIVE « 0400000 RPM/SEC 

AFTERBURNER FUEL FLOW a 940000 


STA 6 


324260 
116148' 
0449424E-01 
194479 
1 0 .000000 
2 4 4066 
O.OOOOOo"""** 
0445519 " 


STA 7 


304377 
2660? 78* 
%5695?7" 
199.576" 


0478756 
*440136 * 


VARIABLE GEOMETRY — 

FVOP « -8.50040 
CVOP a 4 40000 
THROAT AREA a 666400 


(d) Operating point 4 (wet design, point), 
figure 14. - Continued. 


PSHIPT *-.7314866-0* 
CSHIFT *0 !S09492E-04 






*w* 






INPUT BATA 

OPPRATIND POINT NUflBEO 


ORIGINAL PACK IS! 
OF POOR QUALITY 


T*M0 « 


o.oooo seconds 


^pressure 

TWCRATURn 

JLRIVAnJr 

_N^^Ft*Ow' 

J„mun?v: 
stored mass* 
DGRlVAWtT" 
GMLROY DER?” 

dcita"h~*'* 


STA 9 


H . 69 A 0 

518470** 


193,000 


l0 ^SESt° SPOOL c 9200,00 opM 
DERIVATIVE o 0, 160739B-01 RPM/SEC 

MAIN COMBUSTOR FUEL Flow e l .;0000 

DUEO MASS FI ONS- 
LOW PRESSURE e 1.63999 
PRESSURE « llfmo 
OVERBOARD -0.260009 

COHVEROEO STEADY STATE POINT 

TIME ® 0.0000 SECONDS 


STA 13 

»TA 2,8 

J6J0D0 
696,000"**""* 
"t 667691 E-02** 

jjjjqqo 

696, 000 """"" 

J6,5000*""~" 

0J9062Dt>Sr 

107,000 

6.73723"**" 

t 

o.josnin-oi 

1 

| 

“.AA9973E-01 

1 




STA 3 


R67.000 

malJr 

; : 333a«r 
88, Om" 


STA A 


256,000 

oj«ojir"” n 

0 i 660266 

^sosim-os’ 

0J67162 ' 

167,000 


STA A, I 


70.0000 
178 0,00**^ 

oas7mn-oi 

107,000*""" 
J. 66768 

%6l03S2E-0r 

ojimp-or 

73,9001***""" 


HI0H DFRfu?r?K2 0L S H900.0 RPM 
DERIVATIVE - 0.6986S7E~02RP,V3EC 

AFTERBURNER FUEL FLOW a 2.80000 
VARIABLE GEOMETRY — 

FVGP c -2,30000 
CVOP a A, 00000 
THROAT AREA a 860.000 


STA 6 


31.8000 

** ■»«»»* 

1160 , 00 

*466589"" ” 

***’ 

196,960 

«W«. 

' • 3128001 ^ 0 ) 

*** 

4220700-03* 
.322003 


STA 7 


JO. 6000 

jm:;r 

jimT 

ijirjir 

1.06272*’ 

ijAsjr* 

.imir 


rauJEI ^*.23810SE-02 
CSHIFT C 0,211B76E“07 


^PRESSURE 
JEMPERATURE 
^DERIVATIVE 
JASSW' 
^DERIVATIVE 

JI25P ^ 

DERIVATtii^*' 
ENERGY~DER?~ 
DELTA H **** 


STA 2 


1A.6960 

518470 * 

193473** 


STA 13 


36.6585 
696 .650~"~ 
0 : 323A97E“02 

04765630-01* 

J4638I 

;.1678A7E-03* 

0.2188076-01* 


SPEED SPOOL a 9177,87 »PM 
DERIVATIVE « 0.18U75E-0I RPM/SEC 

MAIN COMBUSTOR FUEL UOW a 1>7 0000 
BLEED MASS ElOWS— 

i?Su P 5IS su, 'e « l.AAJOB 
HIGH PRESSURE a 1? ?u« 
OVERBOARD *0.260$68 8369 


STA 2.2 


364078 
697. 838^ 

mriir 


STA 3 

STA 6 

STA 64 

STA 6 

STA 7 

267.625 

1323.56 

-461692 

68.3117 

256.602 

251645 

04657UE-01 

90.0111 

70.1201 

32.0069 

30.8156 

1777.29 

-4913080-01 

107.268 

U61.02 

0.893216E-01 

1969.81 

0.2663350-01* 


m.412 
! 0.226609E-01 
2*23961 
0 , 30S176E-06 
0.200027 

197.612 

[0.915611 

0.661673 

2.66566 

-# 3O5176E-03 

0 , 663730E-03 

** 198366E-03 
**470639 
754762 

1 * 05091 

**429706 

0.7S50AJE-08 

0. 6196170-06 
0.258877E-01 


”* Vf V2!I ccu »rU0l » 11889 $ ddm 

DERIVATIVE # -.699273E-01RPM/SEC 
AFTERBURNER FUEL FLOW a 2.80000 
VARIABLE GEOMETRY — 

PVGP a -2.50000 
CVOP a 6.00000 
THROAT AREA a 560.000 

( 0 ) Operating point 5. 
figure H - concluded. 


PSHIPT a-»887099E-OS 
CSHIFT a 0 .S 61726E-06 3 




3S? 


*a*M*fc J v <-w V ...13 


ORIGINAL P/lQf gn 
W POOR 


Correction 

coefficient 

J 

2 

3 

4 
8 
0 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 


Va.up 

1,0000000 

1,0000000 

1.0000000 

1,0000000 

0. 99618240 

1.0000000 
1. 0028696 

1.0000000 
1.0099678 
1.0024300 
1.0010843 
1.0267143 
1. 0045 977 
0 9986819 

1, *'016937 
1,0226727 
tt 99627388 

1.0000000 

1.0089293 


Figure is. • Correction coefficients 
for dry design point. 



Figure 16. » Computational flow diagram of turboshaft engine. 











Figure 1?* • Transient response of a small turbosl 
multaneous steps In fuel flow and load. (Value' 
to design point) 





